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ABSTRACT 


This research project proposed a modified unit commitment that 
schedules connection and disconnection of generating units in response to 
load A modified generation control is also proposed that controls steam 
units under automatic generation control, fast responding diesels, gas 
turbines and hydro units under a feedforward control, and wind turbine array 
output under a closed loop array control. This modified generation control 
and unit commitment require prediction of trend wind power variation one hour 
ahead and the prediction of error in this trend wind power prediction one 
half hour ahead. An improved method for predicting trend wind speed 
variation is developed. Methods for accurately simulating the wind array 
power from a limited number of wind speed prediction records was developed. 
Finally, two methods for predicting the error in the trend wind power 
prediction were developed. This research provides a foundation for besting 
and evaluating the modified unit commitment and generation control that 
was developed to maintain operating reliability at a greatly reduced 
overall production cost for utilities with wind generation capacity. 
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EXECUTIVE SUMMARY 


The purpose of the research is the development of: 

(1) a modified unit commitment; 

(2) a modified generation control; 

(3) a trend wind power predictor required by both the modified unit 
commitment and generation control procedures developed; 

(4) a wind power error predictor required by both the modified unit 
commitment and generation control procedures. 

These four developments permit one to answer the following two questions which 
are to be addressed by the research to be conducted within the {federal tfiods 
fioetax Besesccb floor U&fcU&Q 

(1) what is the magnitude of the capacity credit that can be assigned to 
wind energy produced by large arrays based on methods for setting 
and meeting the load following and operating reserve requirements 
within a utility's unit commitment. This magnitude of the 
capacity credit assigned to wind will determine the breakeven 
point in terms of 30 year level ized cost in $/KWH that wind energy 
technology must achieve to warrant large scale implementation by 
utilities. 

(2) Develop a generation control Strategy that minimizes the impact of 
large rapid changes in wind array generation that is caused by 
"rotor synchronization" [1, A-3. III-20] of all wind turbines in the 
array for large meteorological event wind speed changes. The 
utility's steam turbine generation is slow responding and cannot 
compensate for these large rapid wind generation changes. Attempt- 
ing to force these units to compensate for these large wind generat- 
ion changes would cause cycling in these units that would expend 
significant fuel, increase maintenance, and possibly reduce unit 
reliability and lifetime. A modified generation control is proposed 
in this research based on an hour ahead prediction of wind power 
change. This modified generation control would utilize these slow 
responding large steam turbine units up to a limit imposed by the 
utility; fast responding diesels, tydros and gas turbines that are 
not presently effectively controlled, and wind turbine array con- 
trol of wind power output as a last resort. 

This modified generation control strategy has been developed to allow 
the utility to determine the level of participation its large steam turbine 
units should have in compensating for large wind generation changes. This 
modified generation control was designed so that oi’ick pickup units provide 
the principal compensation for the large wind cha;~je. The modified quarter 
hour updated unit commitment strategy would continually unload these quick 
pickup units and replace them with standby economic, peaking, and regulating 
units. The purpose of maintaining quick pickup generation is spinning 
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reserve is to maintain adequate spinning reserve and load following margins 
to compensate for large drops in wind generation If the wind generation 
increase exceeds the allocated combination of the response capability of 
large steam turbines under automatic generation control and the level of 
quick pickup generation connected and loaded (that could be unloaded and 
disconnected in 15 minutes by the feedforward generation control to compen- 
sate for the wind generation increase) * the closed loop wind turbines array 
control would reduce the wind generation rate of change to the level the AGC 
and feedforward control could handle. For wind generation decreases* the 
AGC would again compensate for the wind generation change up to the 
capability allocated to tracking wind generation change. The feedforward 
control of quick pickup units would then be capable of connecting and fully 
loading all quick pickup generation within the spinning reserve within 
fifteen minutes to compensate for wind generation decrease. The quarter 
hour unit commitment would schedule their connection and the feedforward 
generation control would set the gain on their governor controls so that 
they would be properly loaded. These quick pickup units* once oonnected* 
would be controlled utilizing the area control error signal used for regu- 
lating units under AGC. The participation factor on the quick pickup units 
would be adjusted by the feedforward generation control to obtain the de- 
sired generation change out of these quick pickup units and thereby prevent 
units under AGC from exceeding the allocated response capability assigned to 
compensating for wind generation. If the predicted wind generation decrease 
is greater than the combined response capability of the feedforward genera- 
tion control and the allocated response capability of the units under AGC* 
then the closed loop array control attempts to build up a back off reserve 
on the wind turbines by clipping wind generation below the level that would 
otherwise be produced given the present wind speed at the particular time. 

This back off reserve is utilized as a cushion so that when the predicted 
drop in wind generation occurs it is not larger than could be handled by AGC 
and feedforward control within 15 minutes. The development of the back off 
reserve is possible due to the hour ahead prediction interval and the fact • 
that the level of wind generation to be clipped is based on the maxiouQ 
predicted wind generation decrease (trend wind power minus the error in the 
trend wind power prediction) and not just the trend wind power change. 

Although the structure of the modified generation control has been 
developed* it has not been evaluated or bested via simulation. The simulation 
of this generation control utilizing predicted trend wind power changes and 
predicted wind powr errors of the typical utilities that could expect to have 
large wind penetrations would be a subject for future research. 

A trend wind power predictor was developed in this research project* The 
research in this project showed: 

(1) time filtering wind speeds caused significant distortion of the 
maximum* minimum* and average values in viod speeds pcedicibiOQ and 
could introduce significant delays; 

(2) time filtering is not required to determine meteorological event 
propagation direction* the reference groups used to predict wind 
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speed at prediction sites in the wind array or propagation delays 
between referenced and prediction sites. This is a change from the 
wind speed prediction method developed in [20]; 

(3) the reference measurement sites should encircle the wind turbine 
duster at a distance of at least 100 miles away from all wind 
turbine dusters. Meteorological events can propagate at speeds 
between 0-100 mph and thus a 100 mile separation allows one or more 
hour ahead trend wind power prediction; 

(4) the reference grouse should not contain storm cell induced cyclic 
variation because such variation is site specific and time varying, 
□sing reference sites with cydic storm induced variation prevents 
prediction of the trend changes in wind speed that are associated 
with the storm front and can be predicted; 

(5) the reference groups used for prediction should change when the wind 
shift associated with an incoming front first affects a particular 
duster of wind turbines. The reference group should change from 
one that is in front of the wind array in the direction of 
propagation of the initial meteorological event to reference sites 
that are in front of the wind turbine array in the propagation 
direction of the incoming event; 

(6) the use of several wind speed reference sites introduces a spatial 
filtering of wind speed variation associated with a meteorological 
event. This spatial filtering associated with the wind speed 
prediction is shown to cause the predicted wind power variation to 
exceed, the actual wind power poduoed by the array by as much as 10- 
20 %; 

(7) several wind prediction sites are required to produce accurate wind 
array power estimates. The error utilizing a single wind speed 
prediction site to simulate a 90 wind turbine array could be as 
large as 100% depending on the prediction site selected within that 
wind turbine array. The error could be reduced to 25% if three 
reference sites are used. The larger the number of pediction sites 
the smaller will be the effect of site specific effects and wind 
speed pediction errors of any pediction site. If the wind speed 
at each wind turbine is not pedicted due to the computational 
burden# then one should select prediction sites so that each 
pediction site is geographically closest to an equal number of wind 
turbine sites. This method of siting wind prediction sites 
minimizes the site specific effects and error of any one wind 
pediction site on the total array power pediction; 

(8) the study of five different methods of simulating wind array power 
variations indicates that there can be significant differences 
between the results obtained using different methods. These 
differences are minimized as the number of wind pediction sites 
increases. No one method of simulating wind array power variation 
will be most accurate for all wind conditions since the magnitude of 
the error and site specific variation at a wind prediction site will 
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vary with the wind conditions. Since each simulation method mini- 
mizes effects of error at specific sites and accentuate error at 
other sites, no one simulation method can give the most accurate 
estimate of true wind array power variation for all wind conditions; 

(9) the magnitude of the wind array power prediction error depends on 
the magnitude of the storm induced cyclic variation and turbulence 
induced wind power variation that can not be predicted using the 
trend wind power predictor. Die error in the trend wind power 
predictor due to the spatial filtering in the wind speed predictor 
also contributes to wind array power prediction error. This wind 
power prediction error can be the magnitude of the capacity of the 
wind array during storms since the large cyclic variations can cause 
cycling between zero and rated array capacity. The error is so large 
because the cyclic variation which can not be predicted using this 
methodology. The wind power prediction error can be kept below 
10% - 25% for other wind conditions if a sufficient number of wind 
power prediction sites are used to simulate the wind array power 
variation. 

A wind power prediction error predictor was also developed. The wind 
speed prediction error was shown to be a zero mean and normal at sites where 
wind speed prediction is successfully accomplished. The wind power 
prediction error was shown to be slowly time varying. Thus, a wind power 
error predictor was proposed that averages the absolute error between the 
actual array power output and the predicted array power over a 15 or 30 
minute period and uses this error estimate to predict power 30 minutes 
aheacL This predicted error is not allowed to be less than 10% of the wind 
array power output since even though the error may become very small for a 
period of time it does reflect the error that can be expected to occur at 
some time in the future. This wind power prediction error predictor was 
thoroughly tested. The error band around the predicted array power was 
shown to effectively band the actual wind power variation. 

The following accomplishments of this research project are unique: 

(1) the development of wind speed prediction for meteorological events 
and turbulence induced variation. Prediction of wind speed based on 
turbulence alone was performed in [2] , but the magnitude of 
turbulence induced variation is so small compared to meteorological 
event variation that it does not require prediction to assure power 
system reliability and econony; 

(2) the development of a method for predicting the error in the wind 
power predictor; 

(3) the development of wind power prediction methods. The assessment of 
different wind power prediction methods, the effect of increasing 
the number of wind power prediction sites in the array, the proper 
siting of these prediction sites and the proper simulation method 
for producing array power variation from severed, wind speed 
prediciton sites were all investigated in the research; 
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(4) investigation of a modified unit commitment procedure that would 
greatly increase the capacity credit given to wind generation. 
Without trend wind power prediction and wind power error prediction, 
a utility would not be able to oonnect or disconnect nonwind 
generation in proportion to predicted wind generation increase or 
decrease respectively. Thus, although one could achieve a capacity 
credit based on I£LF calculations, the operation of the utility 
effectively prevented wind generation from serving any load since no 
nonwind generation capacity is displaced by wind generation. The 
modified unit commitment procedure proposed would increase load 
following and spinning reserve proportional to the magnitude of the 
wind power prediction error. The magnitude of the spinning reserve 
increase at arty time, which is proportional to wind power prediction 
error at that time, is the amount of the wind generation that is not 
allowed to be counted at meeting load due to the lack of perfect 
prediction of wind power variation. Wind power prediction thus 
permits one to provide capacity credit for wind and improvements in 
wind power prediction accuracy increase the capacity credit given to 
wind array power variation The research performed in this project 
is the only published research on modified unit commitment methods 
that can utilize wind power prediction to modify the 24 hour unit 
ccmnitment based on predicted wind generation changes; 

(5) the development of a generation control strategy based on the one 
developed in [8] but that utilizes both the trend wind power 
predictor and the wind power prediction error predictor for both 
meteorological event and turbulence induced variation. The 
generation control utilizes the control philosophy in the priority 
use of automatic generation control feedforward control, and array 
control but incorporates the effects of predicting meteorological 
events and the effects of wind array power prediction error. The 
generation control strategy proposed would satisfy utility 
reliability requirements while simultaneously assuring economic 
operation. Furthermore, the proposed generation control would limit 
the cy cl icing on large steam units that would increase fuel costs, 
increase forced outages, and possibly reduce unit lifetime. 

However, the methodology has not been integrated into an individual package. 
Thus, the capabilities and performances of the modified unit commitment and 
modified generation control can not be fully quantified and validate, and 
shound be evaluated in a future research project. 
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SECTION 1 


INTRODUCTION 

Two fundamental questions that are to be addressed in the Federal tfiod 

Goergy Five Fear Besearcb Flao* 1&35-1&2Q [1] have been investigated and 

partially answered in this research project. These two questions are: 

1. What is the magnitude of capacity credit that can be assigned to wind 
energy produced from large arrays based cn development of methods for 
setting load following and operating reserve levels [1, pg. A4J. The 

methods for setting and meeting operating and load following requirements 
must meet the utility operating reliability standards [19] but have 
tremendous impact on the economic breakeven price of wind energy where 
utilities would likely begin large scale implementation of wind 
development. If there are no capacity credits given to wind because of 
operating reserve and load following requirements, the long term economic 
breakeven point in 30 year level ized cost for wind energy would be 3 /kWh 
[1, A3]. If the load following and operating reserve requirement give 
capacity credit to wind energy, then wind need not be justified solely 
based on fuel displacement and the economic breakeven price of wind 
energy would decrease substantially. This research project develops a 
modified unit commitment procedure based on an hour ahead wind power 
prediction that can provide significant capacity credits that depend on 
the accuracy of the wind array power prediction at any time; 

(2) Develop a geoecatiQO CQDtEQl strategy that minimizes the impact of large 
rapid changes in wind array generation that is caused by "rotor 
synchronization" [1, A-3, III-20] of all wind turbines in the array for 
large meteorological event wind speed changes. The utility's steam 
turbine generation is slow responding and cannot compensate for these 
large rapid wind generation changes. Attempting to force these units to 
compensate for these large wind generation changes would cause cycling in 
these units that would expend significant fuel, increase maintenance, and 
possibly reduce unit reliability and lifetime. A modified generation 
control is proposed in this research based on an hour ahead prediction of 
wind power change. This modified generation control would utilize these 
slow responding large steam turbine units up to a limit composed by the 
utility; fast responding diesels, hydros and gas turbines that are not 
presently effectively controlled, and wind turbine array control of wind 
power output as a last resort. 

There are four major contributions of this research: 

1. Development of a oew modified unit CQ©®it©eot pcoceduce that utilizes an 
hour ahead prediction of trend wind power change and an hour ahead 
prediction of the error in this trend wind power prediction. This new 
unit procedure would provide significant capacity credit for wind energy 
(8% - 30%) based on both the magnitude of the wind power predicted over 
the next hour and the estimate of the error in this prediction; 
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2. Development of a new geoecatiQQ CQOtlQl that minimizes the impact of 
large rapid "rotor synchronized" array wind power changes on luge slow 
responding steam turbine units; utilizes the fast responding diesels, gas 
turbines, and hydro units to provide the primary compensation for these 
large rapid wind energy changes; and utilizes wind array controls to 
reduce wind energy changes only when required to maintain the utility's 
operating reliability; 

3. Development of a viod power ptedictiOO methodology based (a) on 
improvements in the previously developed wind speed prediction 
methodology [20] and (b) development and assessment of alternate methods 
for simulating predicted wind array power from the predicted wind speeds 
at one or more sites in the wind turbine array; 

4. Development of a wind power predictioo error predictor that can estimate 
the magnitude of the error in the prediction of trend wind power 
variation. This wind power prediction error predictor would estimate the 
(a) magnitude of the large cyclic variations in wind array power due to 
passage of storm cells through an array; (b) the magnitude of turbulence 
induced wind power variations in the array; and (c) the magnitude of the 
error in the prediction of trend wind power change. The need to estimate 
the trend wind prediction error is due to the fact that (1) this error is 
very large compared to the error in predicting electric power demand 
(load) over a 24 hour or hour period and (2) this error has major 
impact on the capacity credit assigned to wind energy in setting load 
following and operating reserve requirements in the new unit commitment 
procedure. This error in wind power prediction would also determine 
whether feed forward generation control of gas turbines, diesels, and 
hydro units needs to be utilized over the next 15 minute period and the 
magnitude of the reduction in wind generation change that the closed loop 
control of wind array power should allow. 

The remainder of this section reviews the work performed under these four 
major contributions and where it is presented in this report. 

The new unit commitment procedure, developed and presented in Section 2 
of this report, is a significant extension of a modified unit commitment 
procedure developed by Michigan State thiversity in [3]. The unit commitment 
procedure developed earlier [3] assumed that both front and storm 
meteorological event induced wind power variation could be accurately 
predicted. The results in Section 4 of this report indicate that the trend 
variation in both front and storms can be predicted one hour ahead but that 
the large cyclic variation in storms cannot be predicted using the wind speed 
prediction methods utilized in this research. The unit commitment procedure 
[3] proposed that a minute updated unit commitment procedure could be 
implemented based on an accurate quarter hour ahead prediction of storm 
induced cyclic variation. Moreover, no explicit method was proposed in [3] 
for estimating the error in the hour ahead wind power prediction as is 
developed in Section 2 and 7 of this report. The unit commitment procedure 
proposed in Section 2 would allow the present 24 hour ahead unit commitment to 
scfc'dule connection and disconnection of large steam turbine peaking, and 
quick pickup units based on a 24 hour ahead prediction of load and a 24 hour 
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ahead prediction of diurnal wind power variation Operating, reserve, and 
load following requirements would be set assuming these 24 hour forecasts of 
wind power were accurate and that the wind array power variation due to 
meteorological events was small. A quarter hour updated unit commitment would 
then be utilized to schedule connection of quick pickup and standby economic, 
peaking, and regulating units to compensate for large meteorological event 
jnAi noA w jjxi power variation. This quarter hour updated unit commitment would 
be based on a hour ahead prediction of trend wind power variation as well as 
an estimate of the error in this trend wind power predictor. The error in 
this predictor would be an estimate of the magnitude of the (a) large storm 
induced cyclic power variation, (b) turbulence induced wind power variation, 
and (c) error in the trend wind power prediction. Methods for setting 
spinning reserve, unloadable generation reserve, and load following reserve 
within the quarter hour unit commitment sure developed. A method that requires 
additional research is proposed for setting operating reserve within the 
quarter hour unit commitment. Methods for meeting these reserve requirements 
as well as minim izing production cost, satisfying minimum shutdown and startup 
constraints, and satisfying minimum and maximum generation constraints on 
generators, within this quarter updated hour unit commitment procedure, are 
discussed. 

A new generation oontrol procedure is developed in Section 2 that is a 
significant extension of the one developed by General Electric in [8]. This 
new generation control procedure utilizes the trend wind power prediction as 
well as the estimate of trend wind power prediction error. This new 
generation oontrol procedure permits the utility to decide the maximum load 
following response capability in MW/minute to be devoted to compensating for 
wind power variations. If predicted maximum wind power change exceeds this 
capability, then a feedforward control of fast responding quick pickup units 
(diesels, gas turbines, hydros), are utilized to compensate for these large 
wind power variations. This generation control is coordinated with the 
quarter hour unit commitment so that standby economic, peaking and regulating 
units are connected to replace quick pickup units that are connected and 
loaded by the feedforward generation control to compensate for wind generation 
decreases. The quick pickup generation is unloaded and disconnected if large 
wind generation increases are experienced. These quick pickup units are the 
primary compensation for large wind generation increases because they can be 
connected and loaded or unloaded and disconnected in fifteen minutes. The 
actual control of the level of generation in the feedforward control would 
depend on area control error in a manner similar to that used on regulating 
units on AGC. If the maximum predicted wind power change in fifteen minutes 
exceeds the maximum response rate capability of the combination of quick 
pickup units under feedforward control and the allocated maximum response of 
steam turbine units for wind variation, then closed loop control of array 
power output would limit wind generation change to a level that could be 
handled by automatic generation control and the feedforward generation con- 
trol. Wind generation increase can be limited to any desired value by the 
closed loop wind generation control. The ability to predict trend wind gene- 
ration change one hour ahead and the use of a Msinjuffl possible wind generation 
decrease to determine the rate of generation decrease for the closed loop 
array control allows the development of a backoff reserve that would help 
compensate for large sudden drops in wind speed. 
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The generation control is a very significant extension of the work 
performed in [8] because it utilizes the trend prediction of meteorological 
event wind power variation rather than a crude Davenport spectrum based 
prediction of turbulence induced variation. The use of wind power error 
prediction as well as trend wind power prediction to determine when the 
feedforward control and the closed loop array control is required and to 
determine the amount of the wind array power variation to eliminate via closed 
loop array control are two other significant contributions of the research 
reported in Section 2. Finally, the ability to limit the response of units 
under automatic generation control for wind generation change and the ability 
to coordinate the quarter hour updated unit commitment and feedforward genera- 
tion control to maintain sufficient quick pickup units in spinning reserve and 
response (load following) capability are contributions of this work. 

i 

The third major contribution of this research is the development of a 
wind array power prediction methodology and the assessment of its accuracy and 
limitations. An improved wind speed prediction methodology is described in 
Section 3 of this report along with a review of previous literature on wind 
prediction. The wind speed prediction methodology requires determining the 
directi cxi of propagation of the meteorological event, the speed of propagation 
of the event and thus the delays between reference measurement sites and the 
wind speed prediction sites in the wind turbine cluster. Methods for 
selecting reference groups to insure accurate trend wind speed prediction for 
meteorological events and for changing reference groups for arrival of a front 
are discussed. The need to provide reference measurement sites that encircle 
all wind turbine clusters at a radius of 100 miles is indicated to be required 
to insure hour ahead wind prediction regardless of the direction of 
propagation of the event. Finally, individual prediction site, reference 
group/prediction site, and reference group/ prediction group predictive models 
are described along with the least square procedure for estimating parameters 
of these models. 

The accuracy and limitations of the improved wind speed prediction 
methodology is assessed in Chapter 4. It is shown that the individual site 
predictive model is much more accurate than either the prediction 
site/reference group or prediction group/ reference group models. It is shown 
that filtering the wind speed records seriously distorts the accuracy of the 
wind speed prediction and is not required for determining reference groups, 
direction of propagation of the front, or prediction delays. Our earlier work 
[3] on wind speed prediction utilized filtering and is shown to seriously 
distort the accuracy of the prediction and is not otherwise required to 
enhance the information in the record required to determine propagation 
direction, reference groups, and prediction delays. It is shown that the 
large cyclic variations due to storm cells passing through the prediction 
sites cannot be predicted due to their variation over both time and distance. 
Reference sites that did not contain storm induced variation was shown to much 
more accurately predict the trend change in wind speed at sites regardless of 
whether they experienced storm cell induced cyclic variation or not. It was 
also found that utilizing reference sites that are closer to prediction sites 
proportionately reduces the prediction interval and reduces prediction error. 

A 100 mile separation between the reference sites that encircle the prediction 
sites was seen to be necessary because the speed of propagation of a 
meteorological event can be as high as 100 mph and as low as 0 mph when wind 
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speeds are 0 to 40 mph. 


Hhe development of the wind array power prediction is begun in Section 5. 
lhe use of several reference sites to produce a predicted wind speed record is 
shown to be equivalent to a spatial filtering of the wind speed profile of a 
meteorological event. This spatial filtering is shown to occur by time 
filtering the reference wind speed record used to produce the predicted wind 
speed at sites within the wind array and ultimately the simulated power out of 
a wind turbine array. The actual wind speed measurement records at the 
prediction sites in the array are also filtered and then used to simulate wind 
array power. The filtering of the actual wind speed measurement records 
causes a significant increase of power out of the array when wind speed is 
near rated wind turbine velocity by increasing the average wind speed and 
increasing the period of saturation of the wind turbines. The filtering also 
significantly distorts average minimum and maximum wind power variation. The 
filtering of wind speed before producing the predicted wind speed records that 
are ultimately utilized to produce predicted array power had almost no effect. 
Thus, the spatial filtering caused by utilizing several wind speed 
measurements for wind speed prediction is inferred. A comparison of the 
predicted and actual wind array power variation for unfiltered wind speed 
records indicates that the predicted wind array power is always greater than 
the actual wind array power due to this spatial filtering. The cyclic 
variation in the actual wind power record is not evident in the predicted 
array power record due to the spatial filtering. The error in predicting the 
wind power variation in fronts is 10-20% but can be as large as 100% during 
passage of storm cells since the large cyclic variation cannot be predicted. 

The use of two wind prediction sites to simulate the wind array power is 
also investigated in Section 5 to determine if using more than one prediction 
site to simulate the wind array power will cause any significant difference. 
The large variation in the mean, rms, and shape of the predicted wind speeds 
at the two geographically close wind prediction sites caused rather large 
differences between the array power simulated using a single wind prediction 
site record and an average of both wind prediction site records. 

The development of the wind array power prediction methodology is 
completed in Section 6 by a study of various methods of simulating wind array 
power using multiple wind prediction sites. Method 1 utilizes a single wind 
prediction record. Methods 2 and 3 produce an average wind turbine power 
record by either averaging wind speed before simulating power from a wind 
turbine (Method 2) or by averaging the simulated power out of a wind turbine 
that experiences the predicted wind speed record (Method 3). The wind array 
power record is produced from this average wind turbine power record by sum- 
ming a delayed average wind turbine power record for each wind turbine in the 
array. The delay used is based on the geographical location of that wind 
turbine from the first turbine to be affected by the propagation of the front 
and the speed of propagation of the front. Methods 4, 5, and 6 for simulating 
wind power simulate the power out of subarrays of wind turbines ttat are 
closest to a particular wind speed prediction site. Method 4 utilizes a 
single wind prediction record for a subarray. Methods 5 and 6 utilize the 
closest two wind measurement records to the wind turbines in the subarray and 
produce an average wind turbine power record for the subarray. Method 5 
averages the two wind speed records and then simulates the average wind tur- 
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bine power record for the subarray. Method 6 simulates the power from a wind 
turbine at each wind prediction site and then averages the two wind power 
records for the wind speed prediction sites that are closest to the subarray 
of wind turbines. The power out of a subarray using Methods 4-6 produces the 
total subarray power record by summing a delayed average subarray turbine 
power record for each wind turbine in the subarray using a delay based on the 
distance from the first wind turbine in the subarray to be affected by the 
propagation of the meteorological event. 'Hie total array power record is 
produced by summing the subarray power records. 

The results of simulating wind array power showed that the wind array 
power based on one prediction site could be more than double that produced by 
another wind prediction site in the array, even when the two sites are 
within 5 miles of each other. The wind array power using methods 2-6 varied 
by 25% based on the different number of wind turbines that are affected by 
the predicted wind speed at a particular site. These results suggest that 
one should site wind speed prediction sites so that each prediction site is 
close to an equal number of wind turbines. This poliq- of siting wind 
prediction sites will guarantee that one prediction site, that may 
experience site specific variations which effect an unknown number of 
turbine sites, does not cause significant errors in the power predicted to 
occur from the array. 

The development of methods to estimate wind power prediction error is 
presented in Section 7 of the report. The first method developed assumes the 
error is time invariant and a zero mean normal process. The prediction error 
is plotted for 22 sites on the SESAME array and is shown to have zero mean and 
be time invariant, h statistical test is used to determine whether the error 
is normally distributed. The results indicate that the prediction error is 
normal when the prediction of trend wind speed is accurate. This test was 
performed for both a front and storm front. 

A second method for estimating wind power prediction error averages the 
difference between the actual and predicted wind array power over a 15 or 
thirty minute period and utilizes this average as an estimate of the wind 
power 30 minutes or an hour into the future. This wind power error prediction 
method does not require that the error is zero mean or normal and assumes it 
is slowly time varying. Moreover, the error estimate is never allowed to 
become less than 10% of the power since the fact that the error is temporarily 
snail for a short period does not indicate that the error will be small at 
some future time. The results indicate that the method is accurate in 
estimating an error band around the predicted trend. 

Conclusions of this research project are given in Chapter 8. The unique 
contributions of the research are highlighted. Recommendations for future 
research are also given. 
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SECTION 2 


JUSTIFICATION AND USE OF HIND POWER PREDICTION IN WIT 
CDMUTMEMT AND GENERATION CONTROL 

An analysis and simulation of wind power variations for square and 
rectangular arrays [3,16] was recently made based on wind speed measurements 
and the wind model developed from these wind speed measurements. These 
results indicate the worst case magnitude of wind power change for passage of 
meteorological events could be much larger than any utility could cope with 
and maintain operation. It was shown that the magnitude of wind power changes 
for passage of meteorological events on a single 350 MW array in a 7000 MW 
utility can seriously reduce operating reliability and economy by 
significantly changing the unit commitment, automatic generation control, and 
economic dispatch schedules and operation. Moreover, it was shown that total 
array capacity changes can occur within 10 minutes and can occur repeatedly 
for passage of a front or storm. Finally, it was shown that near total 
capacity power variations can occur simultaneously on different arrays 20-40 
miles apart in the direction of motion of the meteorological event. These 
results clearly indicate that infrequent meteorological events can cause 
serious operating problems on single wind turbine arrays with less than 5% 
penetration. The several 350 MW arrays contemplated in the Pacific Gas and 
Electric and Southern California Edison systems could cause operating problems 
for these utilities for the very infrequent occasion that meteorological 
events occur on these sites. It was also pointed out that there is a need for 
a modified unit commitment and generation control if wind power penetration 
exceeds 5%. Penetrations above 5% appear to be feasible as wind technology 
improves and the installation of large wind turbine arrays increases. A 
discussion of a modified unit commitment and generation control strategy is 
given in Subsection 2.3. 

The effects of turbulence were shown [16] to be quite large on a single 
wind turbine but were shown to cause small variation as a percentage of 
utility capacity for wind turbine array penetrations of 5% on large utilities. 
The difference in the effects of turbulence and meteorological events in terms 
of the magnitude of array wind power variations and thus their effect on a 
utility is due to the fact that: 

(1) the weather map fluctuations associated with energy spectrum below 5 
cycles/hour are generally correlated between sites in an array and have 
relatively larger energy than "gusts." The high correlations make the 
power variations on each wind turbine appear quite similar and thus cause 
large power variations out of the array; 

(2) the turbulence or "gusts" wind speed variation component associated with 
the spectrum above 5 cycles/hour has less energy than the "weather map 
fluctuation” component and is generally uncorrelated between sites. The 
lack of correlation of turbulence between sites generally will cause 
cancellation of wind variation between the different wind turbine sites 
which greatly reduces the turbulence induced power variations out of an 
array. 
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It ia impossible to determine whether there are or are not any 
meteorological events in the energy spectrum of Figure 1 since information is 
lost in the calculation of energy spectrum. However, the conclusion that the 
wind speed variation associated with the energy spectrum below 5 cycles/ hour 
is of concern in operation and control of utilities is valid whether there are 
meteorological events in this spectrum or not. The validity of the concern is 
based on the significant energy of these variations and the high correlation 
between wind turbines in an array for such variation It will be our custom 
to refer to weather map fluctuations as meteorological events in our 
discussion. 

2.1 THE UNIT COMMITMENT AND GENERATION CONTROL PROBLEMS 

Research performed by Michigan State [17,18] and others [9,10,11] has 
shown that large wind power variations from an array of wind turbines can 
cause significant operating problems for a utility. These problems occur 
because a utility's unit commitment and generation control is based on (1) 
handling small cyclic load variation rather than the larger cyclic wind power 
variation and (2) large load trend change can be accurately predicted 24 hours 
ahead but trend wind cyclic wind power variations have not been predicted 
accurately. Utility practice for conventional loads, which can be predicted 
on a 24 hour basis within 2%, has been to cCDOSCt or commit units in advance 
using a "unit commitment" schedule, and to CQOtCQl operating units already 
connected via set point adjustments to unit governors as load varies about the 
predicted value. Since wind power variations exceed load variations on a 
percentage basis, this practice must be modified. Hast cyclic and slow trend 
wind power variations are both large and unpredictable 24 hours ahead of real 
time. Since wind power variations are usually viewed as negative load to the 
utility's unit commitment procedures, which provide fast responding generation 
(load following requirement) and generation reserves (operating reserve), 
there is a unit commitment PCQblem in providing the proper additional reserves 
for wind power variations in unit commitment schedules made 24 hours ahead. 
Trend and cyclic wind power variations due to meteorological events can, for 
wind penetration levels above normal spinning reserve levels (5% of a 
utility's capacity), greatly exceed both the systems spinning reserve, 
unloadable generation and load following capability. This can cause a serious 
reduction in system reliability and a violation of the utility guidelines for 
reliable operation (NERC Minimum Criteria for Operating Reliability) [19]. 

The reduction in operatingg reliability due to large unpredicted wind power 
variation can be alleviated by increasing the load following and spinning 
reserve levels thus reducing or completely eliminating the capacity credit 
allowed to the wind generation capacity. 

Two qqq&qI problems associated with power system operation for wind 
power variations are: 

(1) The utility's automatic generation control will saturate for long periods 
when the total change in wind generation and a simultaneous load change 
will require conventional generation change that exceeds load following 
capability in a ten ninute interval. This problem violates NERC 
performance guidelines [19]. However, it can be eliminated by imposing a 
£atm penetration constraint on the capacity of all r ind turbine 
generators that can be affected by a single thunderstorm front. This 
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farm constraint, which limits wind generation to be less than the spinning 
reserve level in a utility, not only can solve this control problem but 
also would solve the unit commitment problem discussed above, 

(2) Steam turbine uoltft will cycle as a result of simultaneous Iced and 

generation changes that induce frequency deviations that exceed governor 
deadband of conventional units. This continual cycling of units is 
objectionable to generator operators and can cause increased maintenance 
costs, increased forced outage rates and ultimately reduced unit life. 

The cycling of nuclear units is of concern for safety reasons in addition 
to those mentioned above. The cycliiig problem can occur due to a storm 
front sweeping through a wind generator array causing large power 
variations on successive echelons. An echelQQ penetration constraint on 
the capacity of all WTGs in a straight line normal to motion of the 
meteorological event that experience simultaneous wind speed changes will 
eliminate this cycling problem. .The fast cyclic wind variation, which 
lie in a range between 2.7 x 10"* hz and 1.6 x 10" J hz, can be quite 
large and cannot be eliminated by the echelon penetration constraint 
because these cyclic variations come from wind variations in a front or 
storm that affect widely separated echelons or possibly different arrays. 
The cyclic variations, around trend wind speed variations, can be 
compensated by increasing the response capability of automatic generation 
control. 

These two control problems, like the unit commitment problem, result from 
the fact that there are large cyclic and trend wind power variations. The 
difference between the unit commitment problem and the control problem is one 
of providing sufficient generation reserves that can respond rapidly enough in 
the unit commitment and have sufficient control action within the generation 
controls to properly compensate for fast wind power variation. The farm 
penetration constraint acts to limit instantaneous maximum wind power increase 
or decrease so that unit commitment and control can cope with trend and cyclic 
wind variations. 

Tbcee solutions to the unit commitment and generation control problem 
discussed above can now be explained using Figures 2-4. A hypothetical daily 
load curve is used for illustration. The diurral wind generation is shown as 
a constant and the wind generation variation due to meteorological events is 
shown as a set of cyclic and ramp variations. The effective load is shown as 
the difference between the daily load curve and the wind generation. It is 
met using a unit commitment that starts up and shuts down units to provide 
sufficient reserves to insure operating rex lability. The Qpetafciog tfiSfitve is 
composed of both nonspinning reserve and s< inning reserve. ^Piooiog E6SSLW 
is generation connected to the system and running, quick start units such as 
gas or hydro turbines, and all load curtailment capability available to the 
operator. The QQQ&piQQlog teseCYe generation is counted in operating 
reserve but not in spinning reserve. The operating reserve and spinning 
reserve for each of these three solutions to the unit commitment and control 
problem are shown in Figure 2-4. The UQlGQdOble geoecatlOD ESSfcCYe, also 
shown in Figures 2-4, is negative reserve that permits backing off 
conventional generation. Lhits having unloadable generation reserve are 
operated above their minimum generation levels so that wind generation 
increases can be accommodated by reducing conventional generation without 
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Figure 2. Unit coumitment solution that increases operating and 
spinning reserve by the wind turbine array capability. 
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Figure 3. A generation control strategy that shuts down the array 
for passage of meteorological event. 
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Figure 4. A unit commitment and generation control procedure that 
adjusts unit commitment and load following generation 
control capability for wind power variation. 
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tripping units off line. thloadable generation reserves are important to 
preserving operating reliability at times of very low load levels when every 
base loaded unit remains connected if all can be operated above their minimum 
generation loads. Chloadable generation reserve places a constraint that 
indicates these base loaded units must be operated sufficiently above their 
minimum generation levels so that increases in wind generation can be 
accommodated without disconnecting these base loaded units. 

The automatic geoecatioo watEQl (AGO matches the effective generation 
to effective load variation and thus keeps area control error and frequency 
deviations small. The automatic generation control adjusts the generation 
levels of units under the control so that system generation will match total 
system load. The automatic generation control is implemented at a utility 
control center and determines and adjusts the desired generation level of each 
generating unit in the utility. 

The first solution [9] , which adds the capacity of the wind turbine array 
to spinning reserve, unloadable generation and operating reserve, can be 
observed in Figure 2. Spinning reserve and operating reserve on non-wind 
generation unit commitment are maintained at levels that totally ignore the 
presence of the wind generation that reduces the load carried by the units and 
thus increases system spinning and operating reserve* The unloadable 
generation level on non-wind generation unit commitment is modified at night 
when wind generation is available. This allows for wind generation increases 
that equal the total capacity of all wind turbine arrays in the utility. The 
unloadable generation reserve level, shown in Figure 1, is so large at other 
times of the day that the need to accommodate wind generation increases places 
no constraint on the non-wind generation unit commitment. The effective 
generation curve shows that automatic generation control response set without 
consideration of wind generation variation cannot effectively track effective 
load changes during passage of meteorological events. Older these conditions, 
large frequency and area control error deviations occur during passage of 
meteorological events that would continually violate NERC Minimum Criteria for 
Operating Reliability [19]. No adjustment to automatic generation control to 
increase response rate capability was discussed in this solution and so none 
is indicated in Figure 2. The spinning reserve and operating reserve levels 
in this solution to the unit commitment are large. The wind generation is 
added to the spinning and operating reserve levels that would have existed if 
the generation were not present. This solution [9] would give no capacity 
credit to wind generation in meeting load and would require that fuel be 
consumed and operating staffs be maintained for all generating units that 
would have operated if the wind generation was not available. This results 
in commitment of additional units and thus in increased fuel and maintenance 
costs than would be necessary if operating reserve were adjusted in 
accordance with wind generation changes. Additions to wind generation 
capacity must be justified based on displacement of fuel which requires that 
wind technology be capable of producing energy at a 30 year levelized cost 
of 3^/kWh. Allowing capacity credits to wind generation by minimizingg the 
increase in load following and spinning reserve in Solutions 2 and 4, will 
greatly reduce the economic breakeven price for wind technology. 

The second solution [2] would alleviate the large area control error and 
frequency deviations by shotting down the wind turbine arrays during thie 
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passage of meteorological events. The ramp and cyclic variations have been 
shown to be as large as the capacity of all wind turbine arrays, occur in as 
short a period as 10 minutes and possibly cycle with periods of 20 minutes to 
an hour or more, k utility’s automatic generation control must attempt to 
track such variations in order to keep tie lines at prescheduled power flow 
levels. The power flow on tie lines is maintained at prescheduled power flow 
levels to assure that there is sufficient power flow capability to supply the 
inadvertent power for loss of generation contingencies. The solution to 
shutdown the wind turbine arrays during passage of all meteorological events 
does not appear attractive since the wind energy would be lost due to shutdown 
of the array. If accurate prediction of wind power variations for 
meteorological events can be made, the shutdown of the wind turbine arrays and 
the concommitant loss of wind energy can be avoided 

(1) entirely for fronts since the trend wind power variation is accurately 
predicted and the cyclic wind power variation with period of 1-60 minutes 
is not large; 

(2) in part for storms and thunderstorms since the time of arrival and 
departure can be accurately predicted and thus the shutdown of the array 
for long periods before and after the storm nay be avoided. 

The simulation of power from wind turbine arrays that experience severe 
storm and thunderstorm induced wind power variations indicate almost 
simultaneous shutdown of all units in an array due to high speed or low speed 
shutdown logic. The wind turbine units in the array will all startup and very 
quickly reach rated array power levels. This shutdown and startup can be 
experienced several times as storm cells pass through the array. The shutdown 
of the array for the entire period of storm cell activity may be desirable for 
some utilities that would wish to avoid attempting to compensate for the 
repeated shutdown and startup of the array using automatic generation control 
and feedfQtwatd CQQLEQl of fast responding hydro, diesel, and gas turbines. 

This solution [2] of shutting down the arrays for meteorological events 
did not consider the use of wind power prediction to minimize or eliminate the 
need for such shutdown. The shutdown of the wind turbine array was proposed 
as a method of reducing the control problem of compensating for large wind 
power changes as shown in Figure 3. This second solution does not address the 
unit canmitment problem. 

The third solution proposed in [18] addressed both the unit commitment 
and control problem (utilizing no prediction of wind power variation) by 
limiting wind power variation via the farm and echelon penetration constraints 
mentioned earlier. The satisfaction of the farm penetration constraint can be 
observed in smaller levels of wind generation and variation in Figure 4. The 
result is a modification of spinning reserve and unlcadable generation to 
track effective load and to increase these reserves during passage of 
meteorological events as shown in Figure 4. Operating reserve modification 
with wind generation change was not discussed in this study [18] and thus no 
modification from that utilized, when no wind generation is present, is shown 
in Figure 4. thloadable generation may be slightly increased in this solution 
due to wind generation but not equal to the capacity of all wind turbine 
arrays as in the previous solutions discussed. The increase in unlcadable 
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generation would occur if the maximum increase in wind generation exceeded 
maximum first contingency loss of export or load which would violate the 
farm penetration constraint. If the farm penetration constraints were 
satisfied# an increase in unloadable generation would only be needed during 
meteorological events. Spinning reserve would not be increased by the 
capacity of the array as in previous solutions. Spinning reserve would not 
increase significantly if the farm penetration constraint were satisfied. 

The increase in spinning reserve due to wind generation would then only 
be large enough to insure reliable operation for the continual large power 
variations observed for passage of meteorological events. This increase in 
spinning reserve can be observed in Figure 4. Note also the step change in 
spinning reserve lags the step change in wind power output due to a change in 
the wind speed in the array. The step change in spinning reserve is delayed 
from the change in wind speed in order to permit confirmation of apriori wind 
speed forecasts that the change in wind speed will be maintained over the next 
few hours. The reduction in spinning reserve and unloadable generation 
reserve over that in the previous solutions will significantly reduce fuel and 
operating costs. Ihe automatic generation control will not adequately track 
the ramp and cyclic variations due to passage of meteorological events 
although significant reduction in area control error and frequency deviation 
is possible by (1) the addition of load following and spinning reserve 
capability to unit commitment during passage of meteorological events and (2) 
increasing automatic generation control response and response rate capability 
through adjustment of AGC control parameters to exploit these additional 
reserves supplied by unit commitment. Methods for deciding the additions to 
spinning reserve# unloadable generation# and load following reserves for 
continual large wind power variation was very briefly discussed in [18]. No 
detailed procedures were given for determining the magnitude and duration of 
these reserve additions for both unit commitment and for generation control. 

2.2 A MODIFIED UNIT COMMITMENT AND GENERATION CONTROL •STRATEGY 

Utilities have developed unit commitment programs that can determine 
which units should be connected during each hour of predicted load increase 
and which units should be disconnected during periods of predicted load 
decrease. The selection of the units to be connected and generating or 
disconnected as well as the ordering of the units to be connected or 
disconnected within each hour is based on minimizing fuel costs and startup 
and shutdown costs while satisfying constraints on the minimum length of time 
a unit should remain connected or disconnected# minimum and maximum generation 
levels for each unit# etc. 

The unit commitment procedure can accurately and economically schedule 
the connection and disconnection of generating units 24 hours ahead because 
the load demand can be accurately predicted (within 2%) 24 hours ahead. 
Utilities must limit total wind generation capacity to be less than the 
minimum spinning reserve and unloadable generation levels if the wind power is 
not predicted. The spinning reserve and unloadable generation reserve is 
intended to provide the reserves for loss of generation or export 
contingencies and errors in predicting load variation. If wind power 
variation is not predicted 24 hours ahead, one must consider that the wind 
power is not available in scheduling non-wind generation to assure that proper 
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levels of operating reliability are maintained since wind power levels of the 
array are not constant or follow a consistent pattern over every 24 hour 
period, Thus, the procedures in solutions one to three that either add wind 
generation levels to spinning reserve and operating reserve or limit the wind 
generation capacity to be less than the minimum of spinning reserve and 
unloadable generation reserves is required to maintain operating reliability 
if wind power prediction is not undertaken. Increasing spinning reserve and 
operating reserve by the wind generation capacity gives wind generation no 
capacity credit and thus the wind generation would not be counted as serving 
any portion of the customers load. The only economic benefits for installing 
wind would thus be the value of the fuel replacement on conventional units due 
to wind generation, which requires high fuel costs to justify installation of 
wind generation or reduction in the 30 year levelized cost of wind energy to 3^ 
/kWh. The third solution, limiting wind generation capacity to be less than 
the minimum of spinning reserve arid unloadable generation, would decrease 
connected generation in proportion to the magnitude of the wind generation and 
may only slightly increase spinning reserve levels during periods when the 
array is experiencing passage of meteorological events. The operating reserve 
may be decreased by the capacity credit (8%-30%) of wind array capacity 
assigned to the wind generation although no credit was given in Figure 4 or 
discussed in the research [18] where this solution was proposed. This 
third solution has three major disadvantages: 

(1) the wind generation capacity is limited to the load following capability 
provided by the unit commitment and generation control; 

(2) a possible reduction in possible operating reliability by utilizing 
spinning reserve and unloadable generation reserve, which is intended to 
cover loss of generation or export contingencies and error in load 
prediction, to cover the wind power variation out of all wind turbine 
arrays; 

(3) very little (0%-30%) reduction in operating reserve. 

A modified unit commitment procedure that would update the 24 unit 
commitment every quarter hour based on one or more hour ahead prediction and 
24 hour ahead forecast of wind power variation could overcome all of these 
disadvantages and allow large wind generation penetrations. 

The 24 hour ahead unit commitment schedule would be developed based on 
the 24 hour ahead prediction of load and a 24 hour ahead forecast of wind 
power variation. Methods for forecasting wind power [6] 24 hours ahead could 
accurately estimate diurnal wind power variations. The methods developed in 
[6], however, could not accurately estimate meteorological event induced 
changes. 

A major modifies ton of the 24 hour ahead unit commitment was proposed in 
[3] that would update the 24 hour unit commitment utilizing a one or more hour 
ahead prediction of wind power variation. The predicted trend wind power and 
the error in this wind power prediction would be predicted for one or more 
hours ahead of real time. The predicted wind power will be shown to 
accurately capture the trend change in wind power variation due to fronts, 
storms, thunderstorms, or stationary highs. The error in the wind power 
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prediction will be shown in Section 7 to be at least 10% for fronts and 
stationary highs and could be as large as 100% of the wind turbine array 
capacity during storms and thunderstorms. Although the trend change over one 
or more hours in wind power due to the storm front can be predicted, the very 
large cyclic wind power variations due to passage of individual storm or 
thunderstorm cells through wind turbine clusters cannot be predicted one hour 
ahead since their formation changes over time, and movement was shown to be 
erratic. Die large cyclic variation in these storm cells do not appear to 
be highly correlated between sites that experience the same storm cell and 
thus prediction appears to be difficult if not impossible. The estimation 
of an error band around the hour ahead trend wind power prediction does 
appear feasible from results in Section 7 and is essential in order to 
compute the operating reserve, spinning reserve, unloadable generation 
reserve, and load following capability to be provided within the quarter 
hour unit commitment. Since load power prediction error 24 hours ahead is 
far smaller (2%) than the levels of spinning reserve and load following 
capability in the 24 hour unit commitment schedules for utilities without 
wind penetration, the error in load power prediction is often ignored in 
setting spinning reserve levels in a particular utility even though NERC 
guidelines for operating reliability indicates spinning reserve levels can 
be increased based on errors in load prediction. Since the error for wind 
power prediction is at least 10% for fronts and stationary high and can be 
100% for storms or thunderstorms, operating reserve, spinning reserve, 
unloadable generation reserve, and load following capability should all be 
adjusted within the quarter hour updated unit commitment based on the 
changes in the wind power prediction error that are predicted to occur. 

These reserve levels should be adjusted directly proportional to the size of 
this wind power prediction error level, which reflects (1) the bias in 
predicting trend wind power variation over an hour, (2) the effects of 
turbulence (wind power variations with periods of less than 10 minutes), and 
(3) meteorological event (front, stationary high, storm, thunderstorm) in- 
duced cyclic variation with periods of less than an hour. The operating 
reserve and spinning reserve would thus be adjusted based on the predicted 
trend minus the trend error variation since one must be prepared to commit 
non-wind generation capacity based on the lowest predicted level of win] 
generation. The unloadable generation reserve would be adjusted based on 
the predicted trend plus the trend error variation since one must be pre- 
pared to disconnect non-wind generation or dump wind generation during 
periods of light load, when all non- wind generating units are at or near 
their minimum capacity levels, based on the maximum predicted level of wind 
generation. 

The quarter hour updated unit commitment, that utilizes this one or more 
hour ahead prediction of wind power and error in wind power prediction, would 
allow (1) increased wind power penetration, (2) adjustment of spinning 
reserve, unloadable generation, and load following within the unit commitment 
based on both the wind power prediction and its error, and (3) adjustment of 
operating reserve in proportion to predicted wind generation and predicted 
wind generation error as shown in Figure 5. 

The quarter hour updated unit commitment would schedule connection and 
disconnection of peaking, regulating, and economic units that are available to 
be connected or disconnected over the hour prediction interval. The quarter 
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Time in Hours 


Figure 5. A modified unit commitment and generation control procedure 
that utilir^d predicted wind power variation to adjust unit 
commitment and load following generation control capability. 
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hour updated unit commitment would also connect quick pickup units (hydro and 
gas turbine, and diesels) in order to supply the load following capability 
required above that in the units under automatic generation control to cover 
the predicted load minus wind power change over the next hour. The quarter 
hour unit commitment will attempt to keep the quick pickup units in spinning 
reserve by either not allowing them to be heavily loaded if connected or by 
keeping them disconnected in order to cover the unexpected wind power changes 
reflected in the wind power prediction error. The quick pickup wits may be 
connected and loaded in 15 minutes to ewer wind power variations when these 
variations cause the actual wind power to be less than the predicted trend 
wind power. The quarter hour unit commitment would then connect peaking, 
regulating or economic units on standby in order to unload quick pickup wits 
if the constraints on minimum shutdown or operating time on these units are 
satisfied. -The economic peaking, and regulating units may also be connected 
or disconnected in response to predicted trend wind power variation. 
Interruptible load would not be disconnected if quick pickup and the standby 
economic, peaking, and regulating units along with wind turbine array control 
can supply the load following capability needed. Interruptible load will be 
counted as part of spinning reserve along with disconnected or connected but 
unloaded quick pickup units. Inclusion of interruptible load and quick pickup 
in spinning reserve is common in Europe but generally has not been practiced 
in the NERC procedures. Since interruptible load can be disconnected and 
quick pickup can be connected and loaded in 15 minutes, they can be available 
when needed with a 15 minute prediction interval thus satisfying the 
requirements derived of spinning reserve. 

The hour ahead prediction interval for trend and the error in predicting 
this trend wind power variation is required since it takes at least one hour 
to connect the peaking, regulation, and standby economic units that are 
available to the quarter hour updated unit commitment. The quarter hour 
update interval for this quarter hour updated wit commitment is chosen so 
that quick pickup units can be connected with no more than a fifteen minute 
delay. NERC guidelines, that require mismatch in generation and load be 
alleviated in 10 to 15 minutes, could thus be satisfied. The hour ahead 
prediction interval for trend wind power prediction error is also required to 
permit the preparation of peaking, regulation, and standby economic units so 
they could be connected to replace quick pickup wits, that are connected, and 
disconnected in response to actual wind power prediction error. These quick 
pickup wits that are connected would be controlled to compensate for the 
variation in wind power and would be disconnected as the regulating, peaking, 
and economic units prepared besed on the hour atmd trend and trend error 
prediction were connected to replace them and perform the generation control 
functiorv thlouding quick pickup units increases spinning reserve and makes 
the quick pickup units available to compensate for the error in wind power 
prediction in the next hour. 

The solution to the CQQtfol problem proposed in this research should only 
be concerned with the fast trend and cyclic components of wind power variation 
because the slow trend (diurnal) wind component can be accurately predicted 
and handled as the slow trend load component via normal 24 hour unit 
commitment and economic dispatch. The control solution proposes to utilize: 

(a) automatic generation control that would better track the hour ahead 
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predicted trend up to a prespecified limit; 

(b) a supplementary automatic generation control of the peaking, regulating, 
quick pickup units committed by the quarter-hour updated unit commitment 
to respond to the predicted trend and cyclic variation. These units have 
a fast response that either is not utilized fully or is not included in 
present automatic generation control strategies. This is called feedfot- 
wtd geoeratiQQ central in [8]; 

(c) a coordinated blade pitch control on all wind turbines in single or 
multiple arrays that can clip predicted cyclic wind power variation and 
smooth rapid hour ahead predicted trend changes that cannot be easily 
handled by automatic generation oontrol or the feedforward generation 
control. This is called feedback attay central in [6] . 

The coordination of these three controls would be permitted through the 
hour ahead prediction of trend and quarter-hour ahead prediction of cyclic 
wind power variation. The modified generation control has more than ample 
control capability for tracking the very large cyclic and trend variation 
which could be expected when wind penetrations range from 5-15% of a utility's 
capacity. 

2.3 OPERATING RESERVE REQUIREMENTS ON UNIT COMMITMENT 

The purpose of the operating reserve, as stated in its definition in the 
NERC Minimum Criteria for Operating Reliability [19] , is to provide sufficient 
reserve above firm system load to provide for: regulation within the hour to 

cover minute to minute variations (load or generation), load forecasting 
error, loss of equipment (generation or transmissions), and local area 
protection. The operating reserve is set for each hour and the system unit 
commitment is required to provide this required level of operating reserve for 
each hour of the day. The operating reserve can be split into spinning and 
nonspinning reserve components. 

The specific operating reserve level OR^ for each hour must be provided 
by either units presently committed to base, economic dispatch or regulation 
functions or by quick pickup units such as gas turbines, diesels, hydro, 
pumped storage, battery, or other short term storage alternatives. The total 
operating reserve must be distributed throughout the utility so that each 
local area is adequately protected for contingencies such as loss of generation 
or transmission. A constraint [21] on power system economic dispatch and unit 
conmitment at each hour k that incorporates the above considerations is: 

I(k) + QP(k) + EminfDq - P^kbMDR^ * OR(k) 
i 

I(k) the capacity of interruptible load via contract with the customer 

QP(k) capacity of all quick pickup and storage that could be brought on 

line and used within 10-60 minutes at hour k 

OR(k) the operating reserve requirement set based on reliability method 
for hour k 
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the generation level of unit i at hour k 

the desired capacity of unit i 

the maximal operating reserve allowed on unit i 


Pi(k) 

K i 

HDR^ 


The utility has alternatives that could increase operating reserves if 
the above constraints were violated and the system was considered in an 
emergency condition. The values of DC^ and MOR^ could be increased, voltage 
reduction, public appeal via radio and television to limit use, and load 
shedding could all be used to increase the operating reserve. The latter 
three methods of reducing load are of increasing severity and are not utilized 
without care due to their consequences. 

The effects of large wind generation capacity on the prediction of hourly 
operating reserve capacity over a 24-hour period can be quite large. The 
methods [22J used to incorporate wind power generation require prediction of 
wind power averages for each hour of the day based on hourly wind speed 
average measured at a single site. The predicted wind power record is then 
used to determine the load that must be served by conventional generation. 

This load minus the slow trend in wind generation can be used to determine a 
load duration curve and the operating reserve required to maintain LOLP below 
a certain level. 

This procedure is quite satisfactory for setting operating reserve and 
OR(k) for slow trend and diurnal wind power variations and load variations 
that can be predicted 24 hours ahead and are repeated on a daily cycle with 
minor variations. If this procedure was used in solutions 1-3 for the 
commitment problem, the operating reserve would have been reduced in 
proportion to wind generation over the 24 hour period. 

The above procedure could not be used in the quarter-hourly updated unit 
commitment because it would not allow quarter- hourly update of the operating 
reserve level based on the magnitude of both the load and the fast trend and 
cyclic wind power variations due to significant changes in wind speed for 
passage of meteorological events. If the effects of fast trend and cyclic 
variations were included in this method, the operating reserve levels 
calculated would effectively ignore wind generation and add the wind turbine 
array capacity to operating reserve [9]. This is exactly what was indicated 
in the first three solutions proposed for the unit commitment problem in 
Section 1. 

It should be noted that the wind power record from a single wind turbine 
used in [9] was much more oscillatory than out of an array due to large 
effects of turbulence. The large turbulence induced variation on a single 
wind turbine are uncorrelated between wind turbines and thus are averaged out 
in the total power from an array and are small [3]. The large oscillations in 
the wind power record used in [9] could be viewed as occurring from 
meteorological events based on results in [3] The results in [9] then 
indicate that the wind power generation can be ignored if the above [22] 
method is used to set operating reserves for wind variations that include 
passage of meteorological events. The operating reserve was set ignoring wind 
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power generation in Figures 2-4. However, the meteorological events shown in 
Figures 2-4 occur only twice; during the late afternoon peak in load and 
during the night time low. There is an average wind generation level that 
would reduce OR(k) in the 24 hour unit commitment using the procedure [22] 
described above. Thus, the effects of average wind power on OR(k) in the 24 
hour unit conmitment do not appear in Figures 2-4. 

The above procedure may at times reduce OR(k) above load L(t) from the 
values shown in Figures 2-4 but would allow operating reserve to instantly and 
proportionately increase with increasing wind generation and vice versa. 
Operating reserve QR(k + j/4) in the quarter-hourly unit commitment should 
change with significant changes in wind speed or passage of meteorological 
events as shown in Figure 5 to reflect changes in effective load L(t) - W(t) 
and the variability of wind prediction W(t) for passage of meteorological 
events. No method exists at present that can properly determine proportional 
changes in operating reserve level OR(k + j/4) in the quarter-hourly update 
unit commitment for significant changes in wind speed or passage of 
meteorological events. The transient operating reserve model [24] appears to 
be a promising approach to properly adjusting operating rejerve levels if a 
markov state model could be determined for trend and cyclic wind variations 
that utilize statistics obtained from the 1 hour ahead prediction of trend and 
cyclic wind power variations. 

2.4 SPINNING RESERVE, INAADABLE GENERATION, AND LOAD FOLLOWING REQUIREMENTS 

IN THE fCDIFIED UNIT COMMITMENT 

Spinning reserve, unlaadable generation reserve, and load following 
reserve requirements are discussed in this section. 

The spinning reserve, unloadable generation, and load following 
requirements include a trend and a cyclic component based on the hour ahead 
prediction of trend and trend error of wind power variation. The spinning 
reserve SR(k), unloadable generation UG(k) and load following reserve require- 
ments are set based on the following formulas: 


SR(k) « max {Dr (k) + (L^ - L k )T + - (W k+1 - W k )T + Q^; 0} (1) 

UG(k) « max{D c (k) - (L^ - L k )T + (£ k + (w k+1 - W k )T + 0} (2) 

LF(k) « max{UG(k), SR(k)} (3) 

D R (k),(D c (k)) the maximum first contingency loss of reserve (commitment) or 
increase (decrease) in wind generation at hour k (megawatts) 


L k the 24-hour ahead predicted load at hour k (megawatts) 

W k+ i the hour ahead predicted trend wind generation at hour k+1 made 

at hour k (megawatts) 


T 


.1667 hours/bours - fraction of an hour 


- V - «Vn - w k> T 
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tbe predicted effective Iced change in ten minutes during 
(k, k+1) (megawatts) 

tbe effects of load forecasting error and minute by minute load 
variation above (below) trend load variation that requires 
regulation (megawatt a) 

tbe effect of trend wind power forecasting error, turbulence, 
and meteorological events below (above) tbe predicted trend 
<Vl - W k )T + (megawatts) 

The unit commitment can meet these spinning reserve, unloadable generation, 
and load following requirements through components from each generator 
connected, quick pickup units and interruptible load. The constraints on unit 
commitment for spinning reserve, unloadable generation, and load following 
are: 

a^Kk) + 6 k QP(k) + I min{Dq - P^k), MSq} > SR(k) (4) 

is A 

I(k) the capacity of interruptible load via contract with tbe 

customer at hour k (megawatts) 

QP(k) the capacity of all quick pickup and storage that could be 

brought on line in 10-60 minutes (megawatts) 

A set of generators connected to the trananission grid 

a L percentage of interruptible load counted in spinning reserve at 

hour k 

3 It percentage of quick pickup capacity in operating reserve 

counted in spinning reserve at hour k 

Dq desired maximum generation level of unit i (megawatts) 

P^(k) generation level of generator i at hour k (megawatts) 

MSq maximum spinning reserve level allowed on unit i (megawatts) 

ql(k) + y k QP(k) + i min(P i (k) - Mq,- Muq) > UG(k) (5) 

ieA 

percentage of interruptible load actually interrupted at hour k 

percentage of quick pickup capacity that could be unloaded at 
hour k 

Mq minimum desired generation level on generator i (megawatts) 

MUq maximum unloadable generation allowed on unit i (megawatts) 
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( 6 ) 


lQRjPj^k) > LF(k) 
ieA 


rate of response in M4/ain of generator 1 

Note that these oonstraints allow use of quick pickup and interruptible 
lead to be counted in spinning reserve and unloadable generation as required 
in quartern-hourly updated unit commitments The expressions for setting 
(1,2,3) aid meeting (4,5,6} spinning reserve, unloadable generation and load 
following are based on hourly updates (k) because such updates are those for 
the normal 24 hour unit commitment. Hie variables such as Wu, QkA etc. in 
(1,2,3) must be specified every hour. These same expressions ( 1 , 2 , 3 , 4 , 5 , 6) 
will be used for setting and meeting spinning reserve, unloadable generation 
and load following requirements in the quarter hour (k + j/4; j ■ 0,1 ,2,3) 
updated unit commitment. 


2.5 COMPUTATION, JUSTIFICATION, AM) UTOATE OF SPINNING RESERVE, UNLOADABLE 
GENERATION, AM) LOAD FOLLOWING REQUIREMENTS Q DARTER-HOUR UNIT COMMITMENTS 


The quarter-hourly updated unit commitment requires setting W^r W^, 

Qui., qL at quarter-hourly intervals in or^'er to set SR(k + j/4 - 1) , CG(k + 
j/4 - l) and LF(k + j/4 - 1) in (1), (2), and (3) for constraints (4), (5), 
and (6) respectively. The constants ■ W(k + j/4), W k ■ W(k + j/4 - 1), 


fyjt. * Q w (k f j/4 - 1) and ■ Q^Clc + j/4 - 1) since the levels over 
-l), k + 1) must be decided oi 


(k ♦ j/4 

k + j) must be decided'on at"k ♦ j/4 - 1 based on a prediction record of 
wind power variation W^(t) for te(k + j/4 - 1, k + j/4). Spinning reserve, 
unloadable generation, and load following levels are likewise updated at k + 


- 1 to cover the period (k + - 1), k + -"—) for the quarter hour 

4 4 4 

updated unit commitment computed at k 1 for any k ■ 1,2 ,...,24, j 


« 0,1 ,2 ,3. 


The measures W(k + j/4 - 1) , W(k + j/4) , Q w (k + j/4 - 1) and Q™(k + j/4 - 
1) used in setting spinning reserve (1) and unloadable generation (2) can be 
illustrated in part by Figures 6 and 7a. In Figure 6, the load L(t) i~ 
constant and shows no variation. The wind variation is represented by a ramp 
increase over 5 hours leveling off to a constant level. There is qq etc Qt in 
predicting cyclic or trend wind power variation and qq cyclic Wiod QI load 


Y&ciatiQQ SO that $«(k + j/4 - 1) , Q™(k + j/4 - 1) , <£ (k ♦ j/4 - 1) , Qf(k + 
j/4 - 1) are zero. The effective load to be cooed with by conventional steam 
generation is also shown in Figure 6. The basic spinning reserve level is D R' 
where d R' is the maximum first contingency loss of generation reserve 
component that is shown as a dotted line that tracks the variation in L(t). 


The actual adjustment of total connected generation capacity Pfl(t) • L(t) + 
is not continuous but occurs at discrete times, which is indicated by the 
staircase. The spinning reserve formula reflects this staircase effect by 
requiring unit commitment supplied capacity to always provide the basic 
reserve plus the change in L(t) - W(t) over ten minutes; i.e. , 


37 



figure 6. Spinning reserve and unloadable generation requirements for trend wind power variation. 


°R + ** ^c+j/4^ ” ^k-l “ W k+j/4-l^ T 

W(t) = W k+ y 4 . x + (W k+ y 4 - W^y^) (t - [k + j/4 - 1]) 

L(t) » I* + (I* j. - I*) (t - k) 
k + j/4-1 i t i k + j/4 

The unloadable generation is seen as negative reserve in Figure 6. The 
basic reserve level is D c , which is the maximum of the maximum first 
contingency loss of load or export from the utility or the maximum first 
contingency increase in wind generation. The unloadable generation formula 
(2) again reflects the need for the basic unloadable generation reserve minus 
the projected change in effective load in ten minutes, i.e., D c - ((L k - W k+ y 4 - 
(Lj^! - W k+ j/^i) ] T. The unloadable generation changes again change in a M 
staircase reflecting discrete time unit commitment changes that always supply 
at least the unloadable generation in the formula (2) . 

Figure 7a is identical to Figure 6 except that large cyclic wind 
variations are imposed on L(t) - W{t). Note then an additional spinning 

reserve Q^(k + j/4 - 1) and unloadable generation reserve Q^(k + j/4 - 1) are 
required for j = 0,1 ,2 ,3 and for k as long as the cyclic variation persists. 

Note that as the cyclic variation or error in the predicted trend increases, 

the values of Q^(k + j/4 - 1) and Q^(k + j/4 - 1) increase also. 

The methods for predicting trend wind power variation W(t) one hour ahead 
at t * k + j/4 based on measurements at t = k + j/4 - 1 are discussed in 

Sections 5 and 6 of this report. Procedures for predicting the error Q^(k + 

j/4 - 1) and Q^(k + j/4 - 1) were developed in Section 7 of the report. The 
most satisfactory procedures averages the error between the predicted power 
W(t) and the actual wind power %(t) over an interval [k+j/4-1 - T, k + 
j/4 - 1] just before k + j/4 - 1 and utilize this error as a prediction of the 
wind power error at t = k + j/4. A lower limit on this error would prevent 
the predicted error from becoming smaller than this limit. This limit would 
be set based on weather forecasts. This limit would be increased during 
storms or thunderstorms and would be considerably smaller during fronts, 
stationary high, or non-event wind conditions. 

The predicted power W(k + j/4 - 1) lower limit W(k + j/4 - 1) - Q^(k + 

j/4 - 1), and upper limit W(k + j/4 - 1) + Q w (k + j/4 - 1) at t = k + j/4 - 1 
- T based on measurements att = k+j/4-l are plotted in Figures 8a and 8b. 
The results in both figures are based on simulations on wind power at 81 wind 
turbines sited at a density of l/mile^. The wind power was simulated based on 
a predicted wind speed record at a site in the array and based on the actual 
measured wind speed at the site. The predicted trend and predicted error in 
Figure 8a was for a front and the predicted trend and trend error in Figure 8b 
was for a storm front. These results show that the lower limit was 
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considerably smaller for the storm front than for the front and the upper 
limit was higher for the storm front than for the front due to the large 
cyclic variation in storm cells that cannot be predicted by the trend wind 
power predictor. 

The cyclic variation above W(k + j/4 ) shown in Figure 8a and 8b due to 
cyclic variation in fronts and storm fronts respectively can be either 
(dipped by coordinated blade pitch controls on the wind turbines or be 
compensated by the response of quick pickup units such as diesels, gas 
turbines, and hydro units. These quick pickup units would be committed by 
the quarter hour unit commitment in response to this cyclic variation, and 
thus could be used to compensate for it within a feedforward generation 
control. 

The alternative of clipping of the wind power variation above 
W(k + j/4) by wind turbines would lose "free" wind energy as shown in Figure 
7a but would utilize controls on the units that produce the variation to 
eliminate it. The clipping of this cyclic variation above W(k + j/4) would 
reduce or eliminate the need for unloadable generation reserve by reducing 

or eliminating Q^(k + j/4 - 1) within the closed loop or feedback control 
structure. 

There is no ability to clip wind power variations below the trend W(t) 
unless the wind turbine arrays are scheduled to operate below the hour ahead 
predicted trend W(t). This operation of the blade pitch control to dip 

wind power variation to W(k + j/4) - Q™(k + j/4 - 1) as shown in Figure 7b 
does not eliminate spinning reserve or load following responsibility from 
the quarter-hourly updated unit commitment since the "free" wind power 
dipped below W(t) must be provided by the economic, peaking, or regulating 
units and the startup of quick pickup or disconnection of interruptible load 
by the quarter hour updated unit commitment based on the component of spinning 

reserve and load following requirement reflected in Q^(k + j/4 - 1) . 

It should be noted that the option to clip wind power by part or all of 

Q^(k + j/i - 1) below W(t) or not at all is solely the function of generation 
control and has no effect on unit commitment. Since spinning 

reserve is based on W{k + j/4 - 1) - Q^(k + j/4 - 1), the addition to 
spinning reserve due to wind is 

%(t) - W(t) + oJ(k + j/4 - 1) 

The term P^(t) is the actual power produced from the array at time t. The 
generation control composed of automatic generation control of steam turbine 
units, feedforward control of quick pickup units, and feedback control can 
alltv this addition to spinning reserve to reside on the wind turbines if no 
arra^ control is utilized or to reside on steam turbine, quick pickup unit, 
and wind turbines if feedback control of the wind turbine array power is 


42 



5 | 

E % 















undertaken. 

2.6 MODIFIED GENERATION CONTROL 

The automatic generation control has generally performed the regulation 
function of attempting to maintain frequency close to a nominal value and 
matching generation to load change thus nulling area control error which is 
a measure of the mismatch in load and generation. The load following 
requirement on automatic generation control requires that the parameters of 
the AGC are set to command sufficient generation change in 10 minutes to 
null area control error. The automatic generation control load following 
capability is of no value unless the unit commitment has provided through 
constraint (6) the load following capability given by (3). The load 
following capability within the automatic generation control must also meet 
or exceed that given by (3). Results in [8] indicate that insufficient load 
following capability in either unit commitment or automatic generation 
control will cause large excessive area control errors that are sustained 
for periods much longer than 10 minutes which violates the NERC guidelines 
[19] that require (a) the area control error maxima and average values over 
ten minutes to be below a certain threshold based on system size and (b) 
that area control error must pass through zero in every ten minute period. 

It is conceivable that automatic generation control response capability 
could be adjusted to each update of the quarter-hourly unit commitment to 
provide LF(k + j/4 - 1) given by (3) for the unit commitment. However , the 
use of conventional generation to meet these requirements for wind power 
variation due to significant wind speeds changes or passage of 
meteorological events would: 

(1) increase fuel costs on the units committed to providing LF(k + j/4 - 1) 
in (6); 

(2) increase maintenance costs, increase forced outage rates, and reduce 
unit lifetime due to the large continual cycling of these units. 

An improved generation control strategy would utilize: 

(1) normal automatic generation control that without wind variation is 
totally responsible under normal conditions to maintain system 
electrical frequency at 60 hz and regulate total system generation to 
track load variations; 

(2) array controls that would smooth the effects of turbulence; slow trend; 
fast trend variations due to fronts, storms, and thunderstorms; and 
cyclic variations due to turbulence, fronts, storms, and thunderstorms 
from single and multiple arrays. The closed loop array controls would 
utilize a coordinated blade pitch control of all wind turbines in an 
array based on information about the capability of the utility's 
controls to handle these wind power variation components. Such closed 
loop array controls were discussed ir, [8]. These array controls 
utilized fast wind power variation prediction but make no effort to 
predict power variation from meteorological events since the 
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application of the proposed control does not emphasize meteorological 
events; 

(3) a feedforward generation control developed to control the fast 

responding regulation and peaking units committed by the quarter-hour 
unit commitment update and even faster responding quick pickup units 
that would be committed in the minute unit commitment update. At 
present, quick pickup units are not generally utilized in automatic 
generation control and generally not all peaking and regulation units 
are utilized. If peaking, regulating, and quick pickup units are under 
the automatic generation control, their response rate capabilities are 
generally not fully exploited. Ihus, this feedforward generation 
control would control the units committed by quarter hour updated unit 
commitment. This feedforward control would utilize the response 
capabilities of these units to compensate for normal turbulence and the 
fast trend and cyclic wind variation for meteorological events. 

The normal automatic generation control, closed loop array control, and 
feedforward generation control would be coordinated to provide the best 
control performance needed to maintain reliable operation and minimize the 
total cost of regulation of these wind power variation components. 

The procedure developed in [8] for this modified generation control 
based on prediction of wind power variation at least 20 minutes ahead of 
real time. Although the procedure was evaluated for turbulence (normal) 
wind condition and not for meteorological event wind conditions, the 
procedure could be adapted to utilize the trend wind power and trend error 
wind power prediction as follows: 

(1) determine the ramping capability ap g (MW/min, up and down) for the 
conventional generation on regulation. Then, the maximum ten minute 
change in generation is 

ap g = 10 1 ap g 1 

(2) determine the maximum change in load in a ten minute period (AP^ max) 
based on the 24 hour ahead load prediction 

(3) the maximum wind power change in 10 minutes must satisfy constraints 

1 4P = - 4P w 

based on the NERC operating reliability requirement that the area 
control error must be zero at least once in a 10 minute interval than 
the maximum change in wind that can be accommodated by AGC(Ap UfnJ , v ) must 
also satisfy [8] Wmax 

A W * ap G - * 675 Ap L ’ 22 * 5 

max 

based on NERC operating reliability requirejnent that the average area 
control error be less than some limit (Ld) that is proportional to the 
maximum change in load. The second constraint is more binding than the 
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first. The level of generation change out of steam turbine 

units that a utility Swilling to utilize to compensate for wind power 
change must satisfy the above two constraints but is generally much 
smaller than the level determined by the constraint. The actual level 
of &I ia act so that (1) sufficient AGC response rate capability is 
pcoviaeoto compensate for the maximum first contingency loss of 
generation/ (2) excessive fuel is not consumed in ramping large steam 
turbines to compensate for wind generation changer (3) excessive 
maintenance and unit loss of life is not sustained on these large steam 
units for attempting to compensate for wind power variation. The level 
A set by the utility will be used to determine the feedforward and 
dosed loop control in the next step of the generation control 
strategy. 


(4) given a particular penetration of wind turbines, determine the control 
regime knowing present wind generation P w (k + j/4 - 1) predicted trend 

wind generation W(k + j/4) predicted error (Q^(k + j/4 - 1) 

Q^(k + j/4 - 1)) and the maximum increase P^.(k + j/4 - 1) and decrease 

Pp(k + j/4 - 1) in generation in 10 minutes for units under feedforward 

generation control. 


(i) run open loop without feedforward or dosed loop control if 
> [W(Jcfj/ 4hP»(k+j/ 4-1)] -*Q-(k+j/4-l)»P (k+j/ 4-1) 

w “® x if wik+j^j-itotk+yr-i) >o u 


-4ft jta -(k+j/4-l) < |W(k+j/4) -P w (k+ j/4-1) ] -oi(k+ j/4-l)»P rt (k+j/4-l) 

if W(k+ j/4) -fy(k+ j/4-1) < 0 (7) 

This control requires the allocated change in wind generation for 
normal AGC to exceed the predicted trend change (up or down) plus 
the appropriate error 

(ii) run with feedforward control if 

AP W (k+j/ 4-l)+Pp(k+ j/4-1) > P u (k+j/ 4-1) 

max ** 

if W(k+ j/4) -%(k+ j/4-1) l 0 

a% (k+j/ 4-l)+Pp(k+j/ 4-1) > -P d (k+ j/4-1) 

max 

if W(k+ j/4) -%(k+ j/4-1) < 0 

( 8 ) 


This feedforward control would not be utilized if (7) were 
satisfied. If (7) were violated, the response capability of the 
quick pickup units P F (k + j/4 - 1) would be added to that 
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provided for wind by AGCAP W (k+j/4-1). The quick pickup units would 

max 

be adjusted based on a maximum required change of 
AQP(k+j/4-l) - P u (k+j/4-1) -P u (k+*~-l) -1^ (k+j/4-1) +fy (k+^-l) 

if W(k+j/4)-fy(k+ j/4-1) l 0 

j-1 

P d (k+j/4-1) -P d (k+t-1) (k+y 4-1} +1^ (k+---l) 

if W(k+ j/ 4) -fy(k+j/ 4rl) < 0 

where 

[1- 6(k+ j/4-1) ]QP(k+j/4-l) » [l-8(k+^— l)]QP(k+j/4-l)+ AQP(k+^-l) 

4 4 (9) 

Pp (k+j/4-1) - S(k+ j/4-1) QP(k+ j/4-1) 

Pp(k+j/4~l) - [1- 8(k+ j/4-1) ] QP(k+ j/4-1) 

Pp (k+j/4-1) is the portion of QP(k+j/4-l) that is connected and thus 

can be unloaded in 15 minutes at k+j/4-1. Pp (k+j/4-1) is the portion 
of quick pickup capacity that is either disconnected or connected and 
unleaded at k+j/4-1. It is assumed that quick pickup units can be 
connected and loaded to full capacity in 15 minutes or completely 
unloaded and disconnected in 15 minutes. Hie quarter hour updated unit 

commitment will attempt to reduce Pp(k+ j/4-1) and thus keep the 

response capability Pp(k+j/4-l) to wind generation decreases large 
based on the spinning reserve constraints (1,4) on the unit commitment. 

Reducing Pp(k+j/4-l) by replacing quick pickup units by lower cost 
standby economic, peaking, and regulating units also reduces operating 
costs. Hie actual feedforward generation oontrol would be based on the 
area oontrol error in a manner similar to regulating units on AGC. 

(iii) Run with feedforward and closed loop control if 

Afy (k+ j/ 4-1 ) +Pp (k+j/4-1 ) < P u (k+j/4-1) 

max 

if W(k+j/4)-%(k+ j/4-1) * 0 

(k+j/4-1) +P p (k+j/4-1) < -P d (k+j/4-1) 
max 

if W(k+j/4)-%(k+ j/4-1) 1 0 (10) 
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The wind turbine array control can reduce the wind power increase to 

%(k+j/4-l) -A% (k+j/4-1) +Pp (k+j/4-1) 

max 

S P u (k+j/4-l) (11) 

The wind turbine array control can attempt to limit the wind power 
decrease to 


fy(k+j/4-l) »-[£A%(k+j/4-l) Pp (k+j/4-1) ] 

* P d (k+j/4-l) (12) 

where £ can be either 0 or K n (0 i £g£l) depending on the level of 
the predicted power output. The power decrease can only be re-* 
duced if the hour prediction interval is used to start reducing 
wind power early in anticipation of the total change required in 
the next hour. In addition, if wind power error is large, the 
closed loop array may cause the array to operate at less than its 
uncontrolled open loop generation levels so that sudden shutdown 
or large reductions in wind generation from the array can be 
handled by the AGC and feedforward control. 


This generation aontrol strategy allows the utility to set the level of 
generation change (\P W nw J in 10 minutes out of the steam turbine units on 
automatic generation control and economic dispatch it is willing to devote 
to compensating for wind generation change A P Wnv>y (k+j/4-l) can be changed 
based on the level of wind generation change expected or based on an 
operating policy of the utility operators toward compensating for wind power 
changes. If the maximum increase in Hind generation P u (k+j/4-1) or maximum 
predicted wind generation decrease P d (k+j/4-1) in 10 minutes can be handled 
using the automatic generation control capability AP Wnv > y (k+j/4-l) devoted to 
tracking wind generation, no change in feedforward or closed loop array 
control is needed. However, if the maximum predicted change in wind exceeds 
A P W lwty (k+j/4-l) , feedforward control of quick pickup units is required up to 
a maximum decrease (Pp(k+j/4-l) or increase Pp( k+j/4-1) that is to be 
utilized to compensate for wind generation change. Pp (k+j/4-1) is the 
generation capacity of quick pickup units that are not disconnected or 
connected and unloaded that is counted in spinning reserve and can be 
connected and loaded in 15 minutes. Pp(k+j/4-l) is the quick pickup 
generation connected that can be unloaded in 15 minutes and counted in 
unloadable generation reserve. A QP(k+j/4-l) is the change in the generation 
capacity of quick pickup units connected and depends on the change in wind 
generation minus the change in response of steam units under automatic 
generation control that is affected by connection of standby peaking, 
economic, and regulating units. 


If the maximum wind generation increase P u (k+j/4-1) and decrease 
P d (k+j/4-1) is greater than the combined allocated response of steam unit 
under AGC and quick pickup units under feedforward control, then closed loop 
array control is utilized to reduce worst case wind generation increase to 
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generation incraaae P u (k+j/4-l) includes the predicted trend increase over 
ten minutes plus the error Q^(k+j/4-l) . 

Wind generation decrease must be compensated by nonmind generation 
increase. The additional spinning reserve for wind generation prediction 
error Qu(k*j/4-l) is carried by steam turbine units under automatic 
generation control and by quick pickup units under feedforward control. 

Since ax lman poeible wind generation decrease in ten minutes P d (k+j/4-l) is 
based on the total error Qjf(k+j/4-l) and the predicted hour ahead trend 
change 

w k+j/4 ' 

the wind generation control reduces wind generation earlier than actually 
required and more than actually required to produce a backoff reserve. The 
level of this backoff reserve is based on the value of S 0 selected. This 
backoff reserve means the wind turbine array produces less power than it 
actually can given the wind at any particular time and thus when the wind 
speed actually drops the wind power drop out of the array is substantially 
less. 


It is clear that the proposed unit commitment procedure requires the 
development of the one or more hour trend wind predictor and the error wind 
power predictor to compensate for the inherent delays in starting up 
regulation and peaking, quick pickup units, and wind turbines, respectively. 
The proposed control procedure requires prediction to (1) permit the units 
under automatic generation control and the peaking, regulating, and quick 
pickup units under feedforward generation control to anticipate the large 
trend and cyclic wind power variations predicted and thus effectively 
increase their ability to respond; (2) the closed loop array control to 
anticipate and thus reduce total wind power change and rate of change by (a) 
beginning the wind generation increase or decrease before it actually occurs 
and (b) by clipping cyclic wind power variation making the wind generation 
change easier to cope with by PCC or feedforward generation control; and (3} 
to properly coordinate the system AGC, feedforward generation control, and 
closed loop array control portion of the control task. It should be noted 
that if the array control is capable of anticipating a wind generation 
increase or clipping cyclic wind power variations below the trend wind power 
variation W(k+j/4-l) requires the array to operate below the level possible 
with the wind speeds observed at all wind turbines in the array. If the 
closed loop array controls operate the array below the predicted trend 
W(k+j/4nl), the closed loop array control can compensate for positive or 
s m all negative errors in predicting W(k+j/4-l) since small negative errors 
and all positive errors in predicting wind array power are eliminated. This 
type of control of wind array power will be necessary when errors in 
predicting wind power trend W(t) are large and can change rapidly such as 
during stou« and thunderstorms. 

If either the hour trend and trend error wind power predictor were not 
feasible, the solution to the unit commitment and control problems proposed 
in [17] would be implemented with the following consequences: 
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(1) the maximum wind penetration would be limited by the farm penetration 
constraint to the maximum first contingency loss of conventional 
generation [17]. If the spinning reserve and lead following capability 
were increased with array capacity, significant fuel, operating, and 
maintenance costs, that would be added would significantly hurt the 
economics of wind generation, would be added) 

(2) addition to spinning reserve and load following requirements on unit 
commitment would be required in proportion to the maximum cyclic power 
variations anticipated for the next 24 hour period. These maximum cyclic 
deviations would be for the cyclic wind variation during the worst front, 
storm, or thunderstorm that can be anticipated for that day whether it 
occurs or not. These additional spinning and load following reserves 
would be included in these requirements (3) for the entire day or a 
significant portion of it since the time of arrival of meteorological 
events could not be predicted accurately 24 hours ahead. This addition 
to spinning reserve, unloadabie generation, and load following could be 
significant and again reduce the economic viability of wind generation; 

(3) a response and response rate capability would be provided in excess of 
that required of the automatic generation control and closed loop array 
ooc trol. These oontrol actions increase operating costs since this 
constant adjustment of generation levels increases fuel costs and 
operating and maintenance costs, and since the use of array control to 
smooth total wind power variation reduces the energy output of the array. 
These oosts are in addition to the above costs that exist purely for 
connecting the additional generation since these regulations oosts are 
attributed to continually changing generation levels and costs for the 
loet energy from arrays required to clip or smooth cyclic and trend wind 
power variations; 

(4) the reduoed control performance by lack of feedforward generation control 
units and the lack of anticipation and coordination in the system 
automatic generation control, closed loop array control, and feedforward 
generation centred; 

(5) the operating reserve would not be adjusted for arrival of meteorological 
events but will be set 24 hours ahead based on the worst anticipated 
drops in wind generation over the next 24 hour period. 

The recent HECD decision to install an 8% penetration wind array that 
exceeds typical spinning reserve and load following capability points out the 
need for this new solution [3] to the unit commitment and control problems 
because possibly severe reliability or economic penalties can be anticipated 
if the entire array is built and exceeds the farm penetration constraint. It 
is anticipated that other utilities will eventually desire to install higher 
penetrations (% wind capacity) than typical spinning reserve levels (5%) as 
wind technology improves resulting in larger and more efficient wind turbines 
and thus larger wind generation penetrations. Although there are more 
significant penalties for lack of a modified unit commitment when penetration 
levels exceed the farm penetration constraint, the above economic penalties 
are very significant for wind arrays of any size and thus further development 
of this modified unit commitment and generation control is necessary tc: 
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minimize the economic penalties 


SECTION 3 


WIND POWER PREDICTION METHODS 

Wind power prediction of diurnal trend, turbulence, and meteorological 
event trend wind power variation is required for the modified unit commitment 
and for the modified generation control strategy. The methods required for 
predicting each of these wind power components is different and obviously 
their use in the 24 hour, and quarter updated unit commitment, and in the 
automatic generation control, feedforward generation control, and closed loop 
array control are different. The prediction of the weather map (meteorologi- 
cal event) trend wind power change is by far the most important because 

(1) the weather map fluctuations associated with energy spectrum below 5 
cycles/ hour are generally correlated between sites in an array and teve 
relatively larger energy than "gusts". The high correlations in weather 
map fluctuations make the power variations on each wind turbine appear 
quite similar and thus cause large power variations out of the array; 

(2) the turbulence or "gusts" wind speed variation component associated with 
the spectrum above 5 cycles/hour has less energy than the "weather map 
fluctuation" component and is generally uncorrelated between sites. The 
lack of correlation of turbulence between sites generally will cause 
cancellation of wind variation between the different wind turbine sites 
which greatly reduce the turbulence induced power variations out of an 
array. 

It is impossible to discuss whether there are or are not any 
meteorological events in the spectrum of Figure 1. This information is lost 
in the calculation of the energy spectrum. However, the conclusion that the 
wind speed variation associated with the energy spectrum below 5 hz is of 
concern in operation and control of utilities is valid whether there are 
meteorological events in this spectrum or not. The validity of the concern is 
based on the energy of these variations and the high correlation between wind 
turbines in an array for such variation. It will be our custom to refer to 
weather map fluctuations as meteorological events in our discussion. 

Section 3.1 first reviews previous literature on prediction of diurnal, 
meteorological event, and turbulence induced wind power variation. The least 
squares models used for prediction of meteorological events in this research 
is then presented in Section 3.2. The methodology for properly filtering, 
determining the direction of the meteorological event, and determining 
propagation delay for the event is also presented. 

3.1 REVIEW OF WIND POWER PREDICTION METHODS 

The previous literature on prediction of diurnal trend wind power 
variation and the estimation of turbulence induced wind power variation for 
use in the 24 hour, and quarter hour updated unit commitment is now reviewed. 
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The turbulence prediction method developed in [2] would attempt to 
estimate the peak W Q turbulence induced wind power variation defined by 

P{W(t) i W 0 } - .99 

based on a Kaimal spectrum of wind speed, a model of correlation of wind speed 
between wind turbine sites, and a transformation of wind speed to wind power 
variation for the wind turbine models in the particular array. The estimation 
of the peak wind power W. would also depend on the average wind speed measured 
at the wind turbines in the cluster and on the stability of the meteorological 
conditions at the cluster. The estimation procedure would eliminate the need 
to estimate both the actual wind power and the variation around this estimate. 
Moreover, the Kaimal spectrum used in [2] is considered to be more accurate 
than the Davenport spectrum used in [8]. Finally, since the estimation can be 
updated every quarter hour (or minute if necessary), since the error is 
included in the estimate, and since the actual magnitude of turbulence induced 
variation out of arrays [3] even for meteorological events is small, there is 
no need for prediction of turbulence induced wind power variation. Thus, the 
effects of turbulence for normal or meteorological event wind conditions can 
be estimated and then updated as meteorological conditions change. The 
method [2] is utilized and extended for estimating the error in predicting 
trend wind power variation in Section 2 of this report. This error 
prediction would include turbulence but would also include error in 
predicting trend wind power variation and the large cyclic wind variation 
associated with storms and thunderstorms. 

A methodology for subhour wind forecasts was developed in [7]. The 
approach was developed to provide forecasts of trend 10 minutes ahead for the 
modified generation control as well as one to six hours ahead for the modified 
unit commitment strategies. The OEM method, a regression method, and a 
persistence method were selected for evaluation in [5] based on the following 
criteria for a good predictor: 

* techniques should be easily autcraatable 

* ideally techniques should have some physically meaningful basis 

* any predictors used must be available in real time 

* techniques should be applicable to a variety of forecast output formats 
to meet users needs 

* techniques should be applicable for prediction in time frames ranging 
from 10 minutes to a few hours 

* techniques should permit update to be made easily upon demand 

The mean, standard deviation about the mean, the trend, and the standard 
deviation of the trend component were predicted at successive 10 minute time 
steps from 10-60 minutes ahead using a persistence, an autoregressive, and an 
OEM model. The same four variables were also predicted using the persistence, 
autoregressive, and OEM models for successive hour time steps from 1-6 hours 
and for successive half hour time steps from 1/2 to 3 hours. The results are 
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quite preliminary since the research is at an early stage. Hie persistence 
model and OEM were clearly superior to the autoregressive model for predicting 
all four variables and for all prediction intervals. Persistence performs 
nearly as well as OEM for shorter (fewer iterations) prediction intervals 
using each basic time step (either 10 minutes, 30 minutes, or 1 hour). Trend 
forecasts were generally poor using all three methods and improvements could 
be made if there was a method of discriminating whether there would be speed 
change for a site. The large number of cases in the dependent set, where no 
significant wind speed change occurs and the smaller number of cases where 
change occurs in the set of dependent cases, makes the techniques studied 
relatively less effective in predicting large wind speed changes. The 
research performed in this study suggest: 

(1) knowledge of apriori meteorological information about the arrival of 
meteorological events; 

(2) measurements of wind speed and direction at wind measurement sites that 
encircle the wind turbine cluster and that experience the meteorological 
event? 

(3) measurement of pressure and temperature and their changes at the wind 
measurement sites that encircle the turbine array; 

(4) determination of the speed of the meteorological events from wind speed, 
wind speed direction, pressure and temperature measurements, and their 
gradients over time and space; 

would provide the information required to accurately predict meteorological 
event induced wind power changes. 

A method for forecasting trend wind power variation hourly over a 24 hour 
interval is proposed in [6]. This type of prediction would be useful in 
setting operating reserve, spinning reserve, unloadable generation and load 
following reserve in the 24 hour unit commitment. The model first develops a 
static probabilistic transformation that relates hourly average wind power to 
hourly averaged wind speed at a particular wind turbine site for a particular 
wind turbine model (MOO-2, MOD-1, etc.). This static probabilistic 
transformation of a wind turbine was then used in conjunction with 
semiobjective and model output statistics wind forecasts. The perfomance of 
the wind power forecasts was based on properly forecasting whether average 
wind power output for a MOD-2 lies above or below 600, 1200, 1800, or 2400 
kilowatts. The reliability and skill level for these two wind power forecasts 
was encouraging. 

3.2 A WHO SPEED PREDICTION METHODOLOGY FDR METEOROLOGICAL EVENTS 

An effort was made in [3] to establish the feasibility or predicting wind 
speeds at 26 sites in SESAME array of 27 wind measurement sites in a 80 by 80 
mile area in Oklahoma. The wind speed measurements were taken at a height of 
13 feet and at a sampling rate of one per minute. A correlated echelon model 
was used 

W^t) = mi + A i;j (T) [Wj(t - T) - nij) 
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where 




(T)o. 

j 


03) 


m^rmj means of wind speed at sites i and j over time interval (0,NA) 

04 , 0 .: standard deviations of wind speed at sites i and j over the 

J time interval (0,N ) 


delay between the arrival of meteorological event at site i and 
J J j prediction interval 

p ii( T ii) correlation coefficient of wind speed at site i and the wind 

J J speed at site j delayed by T^j 

A 1 minute sampling period for the wind data 


The correlation echelon model assumes the wind speeds at the two sites i and j 
are both stationary processes that can have different mean and variance due to 
surface roughness and site specific effects. The characteristics of the 
meteorological event captured in the wind speed records is assumed to 
propagate from site i to j with the speed and direction of the motion of the 
meteorological event itself. 


The methodology used to determine the model (3) is 

( 1 ) filter each wind record over time interval [ 0 ,N ] using a moving average 
filter 

( 2 ) calculate m^, m^, o^, and a- of the filtered reference wind measurement 
record VMt) ana the wind measurement record VMt) where prediction is 
desired J 

(3) calculate the correlation 
c ij < T > 

Pufr) - 


13 


a. a. 
i 3 


N XjdcAjX^ (kA— t) 

C^(t) = z -- J - ; JL(k) = VMk) - m4 

k=l N ill 


'13 


(4) find the value of and P^jd^j) that maximizes (x) 

(5) predict W^t) using 

p ij {T ij )a i 

W A (t + Ty) = mi - £Wj (t) - ntjj 
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given record of Wj(t). Determine the mean square error for the predictor based 
on error W^(t) -*W^(t) . 

The results obtained from the SESAME data for ten minute moving average 
filtered data showed that accurate estimates were possible at small 
geographical distances from the reference site. The prediction intervals 
were also very small. The estimation errors were much larger at longer J 
geographical distances. Although there was good quality estimation at small 
geographical distances, the prediction intervals were so small that it was 
questionable whether prediction was actually being accomplished. 

The results obtained for filtering data with an hour moving average 
filter were encouraging because reasonable quality estimation was observed for 
sites reasonably distant from the reference site. The delays for some of 
the sites with reasonable quality estimation was 15 minutes ana "thus the 
possibility that prediction could be performed was indicated. 

A correlation echelon model that would utilize wind speed measurements at 
severed reference sites has the form 


Mi 


a. 

l, 


Wilt) - «i + ^ PylTyJiijWjlt - *ij) - WjJ 


(14) 


where 

Mi 

W ± (t) 

mifiUj 

a i'°j 


p^t) 

T u 


Wjlt-Tij) 


number of sites where measurements are taken 

wind speed estimate at site i 

mean wind speed at site i and sites j = 1,2,..., Mi 

standard deviation of the wind speed at site i and sites j = 

1,2,... ,Mi 

normalized cross correlation of wind speed at site i and site j 

delay between the time meteorological event first effects sites 
j and the time it first effects site i. This delay T- • is 
chosen as the value t where normalized cross correlation (x) 
is maximum - 1 


delayed wind speed measurement record at site j 
A similar recursive least squares model is proposed that has the form 


M i 


w^na) = ^ a^OOWjUN - k i;j ]A) + b L 


(N) 


(15) 


Mi 


i 

where {a^j } i and are chosen to minimize JUjj^a^,... ,3^,5^) = 
j— 1 
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N ^ 

z (W^nA) - z a^CnjW^dn - k^jA) - ^(n )] 2 
n-1 j»l J J J 

where t » Na and T^j ■ k^jA. if the processes ar stationary then 
a ij lN) * m 7 p ij^ T ij^o~ “ a ij 


bi(N) ■ m^ - z a ij®j 


A recursive least squares algorithm requires apriori knowledge of = k^A 

Mf 

but allows (a^(N)} and bdN) to be updated at every sampling period t N « 

- 1 j-1 

N A . The method for selecting = k is critical to the performance of the 
predictor and is discussed in ther next section. 


Three different predictive models were tested. The individual site 
predictive model (15) predicts wind speed at each site i in the wind turbine 
cluster using several individual wind speed measurements j each with its own 
The gcoup&itS predictive model predicts wind speeds at each 
site i in the wind turbine cluster based on (a) an average record 
of wind speed in a group of IL reference sites 


delay T,,. 
prediction 


tut) - i- I w,(t> 

"l * J o 

and (b) the average delay between the prediction site i and the group of 
reference sites 




The group/site prediction model has the form 


w i<t) ■ - T i> + b ic 

The accuracy of this model was shown to be inferior to the individual site 
model but using average wind speed records with average delays from severed 
reference groups was found to greatly improve prediction accuracy over that 
for one reference groups The gtOUpdatOUP predictive model averages wind speed 
records at all prediction sites in a wind turbine cluster and all reference 
sites in a reference group and predicts the averaged wind speed record in the 
wind turbine cluster using a delay averaged over all prediction sites in the 
wind turbine cluster and all sites in the reference groups The averaged 
reference wind record for the Nj reference records is 
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W r (t) 



Wj(t) 


An averaged wind speed record for the sites in prediction group k is 


W k (t) 


-- i 

*k ieL, 


W^t) 


and an average group delay 
iel k 


T(k) 


— - l 


“ft jeJ Q iel k 


between reference site jej Q and prediction sites 

T- • 

13 


The single averaged reference group signal W r (t) is delayed by T(k) to produce 
a recursive least squares predictor of the averaged wind speed record for 
group k. The group/group predictor has the form 


W k (t) * a kr w r (t - T(k) ) + b^ 

The implementation of any of the above models requires 

(1) properly filtering the wind speed records at all sites based on the 
propagation speed of the meteorological event being predicted; 

(2) properly determining the direction of propagation of the meteorological 
event which may or may not be identical with wind speed direction at 
individual wind measurement sites; 

(3) properly determining the delay between the reference wind measurement 
sites and those where prediction is being attempted; 


It is clear that although accurate prediction has been accomplished in all 
cases studied^ the procedure would be much more accurate and be able to be 
implemented more successfully on-line if 

(1) meteorological forecasts of the time of arrival and departure as well as 
the speed and direction of motion of the metoeorological event were 
available; 

(2) measurement of pressure and temperature/ and their gradients were 
available. 

The need to (a) very carefully filter the records/ (b) utilize several 
reference wind speed measurements/ (c) utilize several wind speed measurements 
in the geoegraphical region of the wind turbine cluster where prediction is 
desired/ and (d) compute correlation (Pi->( T ij)) and delay (T^) between all 
pairs of wind measurement sites indicated that wind speed and direction 
measurements are very much corrupted by turbulence and site specific 
variations that mask the meteorological event information. Thus, the need for 
additional measurements and forecasts is apparent. 
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Filtering of the measurement records at reference and prediction sites 
using a 2 minute moving average filter to eliminate turbulence can assist in 
accurately determining the direction and speed of propagation of the 
meteorological event. Filtering the records over longer intervals can cause 
serious distortion in the maximum, minimum, and average of the predicted wind 
speed and cause significant delay in the record. Long filter intervals were 
used in an attempt to predict wind power variations for storm fronts in our 
previous research [20]. Results in Section 4 indicate the filtering 
eliminated the very large cyclic variations and was not needed to determine 
propagation delay or propagation direction for the meteorological events as 
was considered necessary in this previous research. 


The direction of propagation of the meteorological event was not always 
identical to the wind speed direction at individual sites or clusters of wind 
measurement sites. A procedure for determining the direction of propagation 
of the meteorological event was determined by calculating the peak correlation 
F ij ^ T ij) f° r all pairs of wind measurement sites 



l N [W^nA)-*^] [WjUn-k^lAj-nij] 

h ” a . a . 

" n»l it 


1 N 

®i * « V nA) 
" n«l 


o. 

l 


1 

9hi 


z (W^nA) — m ^) 2 
n=l 


If sites i^el are the sites of the system ordered in increasing distance in 

a particular direction a and P io jp( T ioV > p ioio (T ioio ) for ail ^ > io 
indicate site j 0 is advanced by T][p-j 0 'from the record i Q . Thus, P^jq(T i 0 j 0 ) 
being larger than p -joio^ T ioio^ implies that the phenomena hits first arid 
then i 0 and not vice versa since delaying j 0 , where the event effects first, 
achieves the larger correlation. If reference sites contain storm or 
thunderstorm cyclic variations, they cannot be used to determine the storm 
fronts propagation direction or delay since the individual storm cells do not 
propagate at the speed or direction of the storm front. Hie wind speed 
direction at a measurement site can be effected by site specific effects and 
not reflect front propagation direction. Hie observation of wind direction at 
several sites and the above procedure can accurately assess front propagation 
direction. 


Selecting the proper delay T ^ Q - 0 between reference site j 0 and 
prediction site i Q is difficult and cannot be based on a single pair of sites. 
Hie smoothing interval must be properly chosen or the delays between reference 
sites and prediction sites may be meaningless since the delays for different 
reference sites j Q geographically close may give very different values of 
delay T io j p to the site i Q where prediction is desired. Hie procedure to 
determ ine J the delay is to 
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(1) find set of several reference sites jeJ f in a small geographical 
region with a very high correlation p jlj2^ T jlj2^ > • 90# 

(2) find a set of several wind speed measurement sites I in a small 
geographical region (where the wind turbine cluster is located) that 
have high correlations p iii2^ T ili2^ > 

(3) if Tj.i are fairly consistent for all iel jej and J > .60 for all 

pairs J id and jej r utilize the set of reference jej vitn delays to 
predict each site iel. J 

This procedure may not determine the delays of the reference sites 
containing storm induced variation. If one tod accurate forecasts of the speed 
of propagation or possibly other information on pressure and temperature 
gradients, one might be able to obtain the storm front propagation delay and 
thus eliminate the need for the procedure altogether. 

lhe reference sites should be at a radius of 100 miles from the set of 
wind turbine clusters to permit hour ahead wind power prediction. The 
propagation speed and direction of the front is often different than the 
wind speed and direction at a wind measurement site. The propagation speed 
of the front can be either much larger or smaller than the wind speed* 

Since front propagation speed can be as high as 100 miles/hour, the 
reference sites must be at least 100 miles from the wind turbine clusters to 
insure hour ahead wind power prediction. If the frontal propagation speed 
is less than 100 raph, then the prediction interval will be longer. If an 
hour ahead prediction is desired for the unit commitment and generation 
control strategies discussed in Section 2, the wind power prediction would 
be provided before it is needed if the propagation delay was greater than 
one hour. 

The reference sites to be used to predict the wind speed at the wind 
turbines in the array should not contain storm or thunderstorm cyclic 
variations since these cyclic variations are not correlated from site to 
site and thus cannot be predicted in the wind turbine clusters from the 
reference measurements. The use of reference sites without storm activity 
allows prediction of the propagation of the front containing the storm 
cells. The cyclic variation in the storm is then predicted or forecasted by 
predicting the error around the wind prediction itself using a procedure 
discussed in Section 7 of this report. 

A procedure for changing the reference measurement sites and associated 
delays for a wind shift is proposed. The reference sites for a wind turbine 
cluster in a large wind turbine array is changed when the wind shift begins 
to affect the sites in that cluster. This procedure requires that the wind 
speed measurements reoord at reference sites in the direction of the 
incoming front be saved and used to determine delay and the predictor's 
parameters when the reference sites and delays are switched to those for the 
incoming front. 
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SECTION 4 


A MODIFIED WIND SPEED PREDICTION METHOD 


4.1 INTRODUCTION 

The wind speed prediction methodology given in the previous section has 
several modifications from that given in our previous work [20]. These 
modifications are based on the work presented in this section. These 
modifications are: 

(1) a moving average filter with long (10 minute - 2 hour) smoothing 
intervals may at time be necessary to increase the correlation and thus 
determine the meteorological event propagation direction. A moving 
average filter with a short 2 minute smoothing interval must be used 
for determining reference groups and prediction groups of wind 
measurement sites, and prediction delays T^j. 

(2) the smoothing interval used for the measurement records for producing a 
predicted wind speed record should never exceed 2 minutes or otherwise 
the average, maximum, minimum, and shape of the cyclic variations in 
that predicted wind record become greatly distorted; 

(3) use of reference group measurement records that contain storms should 
be avoided since the wind speed prediction errors are much more than if 
reference without storm induced variation are used. 

The error in predicting the propagation of a storm will be large since the 
shape of the cyclic wind speed variation change over time as the storm 
propagates. Since storms are local and since there may be several local 
storms, the prediction methodology will attempt to relate the storm in the 
reference group and any prediction site whether the reference group storm 
has propagated to the prediction site or not. Use of a reference group that 
does not contain the storm will thus more accurately predict the wind speed 
at all prediction sites whether they contain the same storm, a different 
storm, or no storm. 

The meteorological event containing the wind shift and a set of storms 
is analyzed in Section 4.2 using a set of maps containing the wind speed and 
direction at the 27 sites in the SESAME array at a particular time. These 
maps are plotted at 20 minute intervals until the storms appear and at 10 
minute intervals thereafter. An analysis of the effects of filtering on the 
accuracy of the wind prediction record is given in Section 43. An analysis 
of the effects of filtering on the estimation of the prediction delays as 
well as an analysis of the selection of reference sites on the accuracy of 
the wind prediction is given in Section 4.4. 
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4.2 PROPAGATION OP A STORM FRONT 

The methodology foe predicting wind speeds developed in the previous 
section is now tested in data from the SESAME array. This array of 27 
meteorological measurement sites are located in an 80 x 00 mile square area 
near Tulsa, Oklahoma* The data utilized in this study was collected over a 
3 month period in the spring of 1979. The wind speed in longitudinal and 
latitudinal directions, the nondirectional wind speed, pressure, 
temperature, and rainfall were all measured at these sites. The latitudinal 
and longitudinal wind speeds were used to determine wind velocity, 
magnitude, and direction of every site in the array. The nondirectional 
wind speed measurements were also retained, but the other meteorological 
data was unfortunately discarded at an earlier stage of the research [3]. 

This pressure and temperature information could have been quite useful in 
determining and confirming the speed and direction of motion of the 
meteorological event as indicated in Section 3. 

The propagation of a storm front containing the propagation of the wind 
characterized by a sudden short duration drop in wind speed and a set of 
large cyclic wind speed variations in a set of local storm cells is now 
analyzed. The storm front occurs over a period from 3:00 - 10:00 ptm. on 
May 2, 1979. The wind measurement sites do not experience the wind 
direction change until 6:40 p.m. as observed on the map of wind speed and 
direction change plotted in Figure 9 for intervals of 20 minutes up until 
7:00 p.m. and then at intervals of 10 minutes to 10:10 pm. The wind shift 
propagates in a southeasterly direction based on both the direction of the 
wind and the propagation of the wind shift line on Figure 9j to 9t. The 
wind shift has affected all the sites by 8:40 and propagates at a speed of 
40 mph* 

The wind shift is accompanied by a set of recurrent storms that 
propagate in a southeasterly direction and then appear to split into two 
parts. One storm remains near sites 8, 16, and 28 and the other propagates 
in a southwesterly direction. This second storm propagates through the 
entire array in approximately three hours. These two storms disappear after 
9:30 p.m. A second set of storms appears near site 2 at 8:50 p.m. and 
propagates in a southwesterly direction in a path similar to the portion of 
the first storm that continued to propagate through the entire array. Sites 
on the western side of the array are totally unaffected by the storm 
propagation. If reference sites used in prediction and prediction sites, 
where wind speed prediction is desired, both contain the large cyclic wind 
speed variations of a storm, the propagation of the storm will be predicted 
in the results to be presented. If either the reference or prediction sites 
do not contain storms in their wind measurement records, the propagation of 
the wind shift is predicted. 

The storm is evidenced by a distortion of the wind speed direction from 
the direction of propagation of the front or wind shift. The distortion can 
appear as a large circular flow of wind around an eye, where wind speed is 
very small. The eye and associated circular flow will then propagate in a 
direction that appears to be related to but not identical to the direction of 
wind shift propagation. The distortion in wind direction can also be seen 
as a circular flow around both sides of an eye, where wind speed is small. 
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Figure 9a. Wind speed and direction maps for a front for the SESAME array 
for. 20 minute perioas from 3:40 - 10:10 p.m. 
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WIND MAP 1600 05/02/79 
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Figure 9b. Wind speed and direction maps for a front for the SESAME array 
for 20 minute periods . rom 3:40 - 10:10 p.m. 
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Figure 9c. Wind speed and direction maps for a front for the SESAME array 
for 20 minute periods fronr3:40 - 10:10 p.m. 
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Figure 9d. Wind speed and direction maps for a front for the SESAME array 
for 20 minute periods from 3:40 - 10:10 p.m. 
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Figure 9e. Wind speed and direction maps for a front for the SESAME array 
for 20 minute periods from 3:40 - 10:10 p.m. 




68 


«»*•; 


WIND 


MAP 1720 05/02/79 

WIND MAP AT 1720 PM 

♦ ♦ ♦ ♦ ♦ 


♦ 



♦ ♦ ♦ ♦ ♦ ♦ ♦ * 

Figure 9f. Wind speed and direction maps for a front for the SESAME array 
for 20 minute periods from 3:40 - 10:10 p.m. 
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Figure 9g 


Wind speed and direction maps for a front for the SESAME array 
for 20 minute periods from '3:40 - 10:10 p.m. 
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Figure 9h. Wind speed and direction maps for a front for the SESAME array 
for 20 minute periods from 3:40 - 30:10 p.m. 
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Figure 9i. Wind speed and direction maps for a front for the SESAME array 
for 20 minute periods from 3:40 - 10:10 p.m. 
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Figure 3j. Wind speed and direction maps for a front for the SESAME array 
for 20 minute periods from 3:40 - 10:10 p.m. 
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Figure 9k. Wind speed and direction maps for a front for the SESAME array 
for 20 minute periods from 3:40 - 10:10 p.m. 
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Figure 91 



. Wind speed and direction maps for a front for the SESAME array 
for 20 minute periods from 3:40 - 10:10 p.m. 
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Figure 9m. Wind speed and direction maps for a front for the SESAME array 
for 20 minute periods from 3:40 - 10:10 p.m. 
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Figure 9n. Wind speed and direction maps for a front for the SESAME array for 
20 minute periods from 3:40 - 10:10 p.m. 
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Figure 9o. 


Wind sneed and direction maps for a front for the SESAME array for 
20 minute periods from 3:40 - 10:10 p.m. 
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Figure 9p. 



Wind speed and direction snaps for a front for the SESAME array for 
20 minute periods from 3:40 - 10:10 p.m. 
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Figure 9q. Wind speed and direction maps for a front for the SESAME array for 
20 minute periods from 3:40 - 10:10 p.m. 
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Figure 9s. Wind speed and direction maps for a front for the SESAME array for 
20 minute periods from 3:40 - 10:10 p.m. 
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Figure 9t. Wind speed and direction maps for a front for the SESAME array for 
20 minute periods from 3:40 - 10:10 p.m. 


+ 


83 


2030 



+ ♦ 


♦ 


♦ + + 


+ 


+ 


Figure 9u. Wind speed and direction maps for a front for the SESAME array for 
20 minute periods from 3:40 - 10:10 p.m. 
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Figure 9v. 
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Wind speed and direction maps for a front for the SESAME array for 
20 minute periods from 3:40 - 10:10 p.m. 
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Figure 9w. Wind speed and direction maps for a front for the SESAME array for 
20 minute periods from 3:40 - 10:10 p.m. 
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Figure 9x. Wind speed and direction maps for a front for the SESAME array for 
20 minute periods from 3:40 - 10:10 p.m. 
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Figure 9y. Wind speed and direction maps for a front for the SESAME array for 
20 minute periods from 3:40 - 10:10 p.m. 
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Figure 9aa. Wind speed and direction maps for a front for the SESAME array for 
20 minute periods from 3:40 - 10:10 p.m. 
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Figure 9bb. Wind speed and direction maps for a f 
20 minute periods from 3:40 - 10:10 p 
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Figure 9cc. Wind speed and direction .naps for a front for the SESAME array for 
20 minute periods from 3:40 - 10:10 p.m. 


*■ 


92 


2200 


OF* POOR QUALITY 


♦ ♦ ♦ 


♦ ♦ 


♦ 


♦ ♦ ♦ 



♦ 


♦ 


♦ 


♦ 


♦ 


♦ 


♦ 


♦ 


+ 


♦ 


<* 


♦ * * 


*■ * * 


Figure 9dd. Wind speed and direction maps for a front for the SESAME array for 
20 minute periods from 3:40 - 10:10 p.m. 
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Figure 9ee. Wind speed and direction maps for a front for the SESAME array for 
20 minute periods from 3:40 - 10:10 p.m. 
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The aet of tbs wind , can help explain the inconsistencies found in 

application of the win' iction methodology to the 3:00 - 10:00 pum. 

record on Nay 2, 1979 shuwn in Figure 9. The direction of propagation of 
the meteorological event could not be predicted with confidence and the set 
of prediction delays were not id ways consistent with the geographical 
distance from the reference sites using the methods developed in [20]. Both 
of theses inconsistencies can be explained based on the analysis of the wind 


The meteorological event propagation direction of the incoming storm 
front is Known to be from northwest to southeast from the wind maps. Dm 
methodology used to determine meteorological event propagation direction, 
discussed in Section 3, should indicate that the correlations in 

the first five columns should be larger than in the first five rows 
given the ordering of sites in Thble 1 is in the known propagation J J 
direction. Elements corresponding to site 9, 11, 13, 17, 19, and 21 had one 
or more elements Pji&u) in the first five rows larger than the elements in 
the first five ooldmns These six sites did not experience the first storm 
that began at 6:30 as shown in the wind maps in Figure 9. Since the long 
smoothing interval emphasizes the propagation of the large and long duration 
wind speed increase associated with the storms, the correlation for sites 
without storms would be small and would not necessarily indicate the 
propagation direction properly as was observed from the correlation in ThfcLe 
L. Ignoring correlations from site 9, 11, 13, 17, 19, and 21, the 
meteorological event propagation direction would be predicted perfectly 
using the methodology developed. 

The prediction delays given in Table 2 are not consistent with 
geographic distance from reference sites 1-5. Group 1 (6,8,26), Group 2 
(10^27,16), Group 6 (22,23 ,24) and Group 11 appear consistent with 
geographical distance and the propagation of the part of the first storm 
that propagates the entire length of the array. Sites (7,9,11) in Group 3 
experienced the second storm but did not experience the first storm. Die 
long delay between sites (7,9,11) in Group 3 and the reference sites appear 
to be an attempt to correlate the first storm at sites 1-5 that first appears 
at 6:30 and the second storm that first appears at 9:30 in sites 1-5 and 
propagates toward sites in Group 3. With storms being local in nature and 
recurrent, it appears to be difficult to track a storm path and correlate 
the reference group that experiences the same storm as a specific wind 
turbine cluster. The prediction of the cyclic variation in storms is 
difficult due to changes in the shape of the large cyclic wind speed 
variations in storms. Thus, efforts to track storms and to accurately 
predict the wind speed variation caused by the storm in a particular wind 
turbine array appears to be nearly impossible in real time with measurements 
available. 

The delays in sites in Group 4 (12,14,15) , Group 5 (13,17) , Group 7 (19) 
and Group 9 (21) are small and these sites do not experience the first or 
second storms. Since the smoothing interval is long and thus emphasizes the 
effects of the large and long duration wind speed increase associated with 
the storm, the short duration drop in wind speed for the wind shift is 
filtered and time shifted. Thus, accurate prediction of the wind shift is 
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Table 1. Peak correlation Matrix P««(T<<) for 2 hour filtered data fro* 3:00 • 10:00 p.M. on Hay 2, 1979, 
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difficult using the data with this long smoothing interval. The results 
from Table 2 indicate that the delays for these sites in groups 4, 5, 7, and 
9 are short compared to the propagation speed of the wind shift (40 mph) and 
their geographical distance from sites 1-5. 

The errors and delays fOr the individual site, group/site, and 
group/group model are given in Table 3. Note that the error for sites in 
Groups 4, 5/ 7, and 9 that had no storm are small. The errors for site in 
group 6 and 11, that contained the storm that appeared in the reference 
group were not large. Sites in group 1 and 2 contained several storms and 
thus the errors were large since only the propagation of the one storm was 
predicted. 

The error: s for the individual site model was smallest since the 
informaton in each of the individual site records can be used for 
prediction. The errors for the group/site model are generally considerably 
larger than for the individual site model since only one averaged record is 
utilized with a single averaged delay, (fee of several reference groups can 
greatly improve the performance of the group/site model. The group/group 
model estimates the average wind speed at a group of measurement sites in a 
wind turbine cluster rather than the individual sites in that cluster. This 
group/group model that estimates the average wind speed in the wind turbine 
cluster based on average wind speed record of a reference group and an 
average delay is by far the worst model. It not only disregards the very 
significant wind speed average, rms and cyclic variations in each record at 
sites in the wind turbine cluster but also ignores this information at 
individual sites in the reference group Methods of utilizing the 
individual predicted wind speed records in a wind turbine cluster to 
accurately simulate wind power out of the wind turbines in the array is 
discussed in Section 6 of this report. 

4.3 THE EFFECTS OF FILTERING ON WIND SPEED PREDICTION 

The long smoothing interval used to determine the meteorological event 
propagation direction, reference groups, and prediction groups for the 
meteorological event should not be used to produce the actual wind speed 
prediction record. Comparison of the actual and predicted wind speed 
records in Figure 10 with a 2 hour smoothing interval and those in Figure 11 
and 12 with 10 minute and 2 minute filtering intervals shows significant 
differences. It is clear that filtering very significantly modifies the 
actual wind record but has substantially less effect on the predicted wind 
records that begin at 270-290 minutes in Figure 10-12. Note that the peak 
on the actual site 22 record att* 310 minutes (8:10 pm.) was reduced from 
64 to less than 50 if the smoothing interval is 10 minutes rather than 2 
minutes. The higher frequency fluctuations were filtered out of both the 
predicted and actual record by the longer smoothing interval. The magnitude 
of the trend change in predicted speed (slope of the average wind speed 
before t = 300 and slope of the average wind speed after t = 300) on the 
predicted wind speed records at site 22 experienced little change with the 
longer filtering interval. Similar conclusions are reached in the analysis 
of sites 23, 24, and 25 in Figures 11 and 12 for 10 and 2 minute smoothing 
intervals. 
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Site 

No. Averaging 

ORIGINAL PAGE IS 
Op POOR QUALITY 

Reference Group - 

v? '•aginc 

3eference~a7Td '^retfict4ofl 
Grouo Averaging 

.maximum 

Site 

Delay 

Error 

(nroh) 

Group/Site 
Delay 
■ minutes' 

irror 
•oh) . 

Group/(iroup 

Delay 

(minutes'. 

trror 

firnh* 

fi 

6 

1 

1.51 

l 

2.69 

2 

2.18 

■ 

0 

U 

3 

0.49 

1 

1.28 

2 

2.13 

[' 

26 

1 

1.22 

1 

2.84 

2 

2.18 

2 

10 

53 

1.00 

35 

2.41 

30 

2.61 

2 

27 

41 

0.67 

25 

1.89 

30 

2.61 

i 2 

16 

51 

1.06 

31 

2.89 

30 

2.61 

3 

7 

208 

0.51 

201 

C.60 

198 

0.62 

3 

9 

231 

0.23 

179 

1.37 

198 

0.62 

3 

11 

227 

0.28 

215 

0.54 

198 

0.62 

4 

12 

24 

0.43 

11 

2.78 

55 

1.07 

A 

•* 

14 

29 

0.64 

10 

* .41 

55 

1 .07 

4 

15 

338 

. 0.11 

145 

0.52 

55 

1 .07 

3 

13 

200 

0.31 

T23 

T .97 

IT? 

T.33 

5 

17 

198 

0.66 

116 

2.96 

119 

1 .98 

6 

22 

275 

0.36 

269 

2.08 

279 

0.42 

6 

23 

290 

0.15 

284 

0.37 

279 

0.42 

6 

24 

290 

0.33 

236 

0.85 

279 

0.42 

7 

19 

22 

0.75 

9 

2.24 

9 

2.24 

8 

2C 

363 

0.14 

346 

0.39 

346 

0.39 

9 

21 

30 

0.47 

16 

0.37 

16 

0.37 

10 

25 

1 

0.63 

1 

1.72 

1 

1.72 

11 

18 

259 

0.36 

242 

0.35 

242 

0.85 


Table 3. Errors and delays for individual site, group’site, and group/group models 
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Figure 11. Actual a1^~7^dTcted~v^ . p u iu: m )limi , ^ ^ m 

DATA USING 19 MINUTE «OVIN* AVERAGE RECORDS AND SITES 1-5 AS REFERENCE^ 
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Comparison of the actual wind speed records in aitas 22-25 ara no 
further than 15 miles apart shows radically different shspad wind spaed 
cyclic variations. Either the cyclic wind speeds variation within a storm 
must be radically different at different points within the storm; or the 
cyclic variations must change radically with time; or bc^h. Thus, it appears 
impossible to predict these very large cyclic storm induced variations at 
any site because (1) the prediction methodology appears to filter them out 
(and only captures the average wind speed increase asso c ia te d with the 
storm) , (2) the large cyclic wind speeds at earn site are so different from 
other sites dose to it, and (3) the cyclic variations at the prediction 
sites are also different than at the reference sites. 

The 2 hour filtering interval totally distorts the actual and the 
predicted wind speed records. The large cyclic variations in the actual 
records completely disappear and the actual and predicted records appear to 
be nearly identical. The relatively short storm induced pulse at site 25 is 
completely eliminated, which explains why the delay chosen is 1 minute 
rather than 270-290 minutes. The longer storm induced variations on sites 
22-24 are indicated as a 2 hour ramp increase in wind speed starting at t ■ 
180 (6:00 pm.) and reaching a maximum between 300 (8:00 pun.) and 330 (8:30 
ptsL) when the actual sharp peak in wind speeds occur. The wind speed at 
maximum (near t - 180) is larger and the average wind speed during the storm 
(t > 300) is smaller losing the 2 hour filtering. 

The errors between the actual and predicted wind for 2 minute, 10 
minute, and 2 hour smoothed data are given in Table 4. Note that the errors 
increase dramatically as the smoothing interval is shortened and the size of 
the error is small at sites that have no storm (13,17,19,21,12,14,15) in 
the 2 minute filtered data. 

4.4 THE EFFECTS OF (SINS DIFFERENT REFERENCE GROUPS 

The prediction of wind speed at sites 22-24 captured the propagation of 
the storm observed in reference group 5. The storm induced cyclic variation 
changed so much over time that the prediction errors were quite large for 
the 2 minute data despite the fact that the storm propagation had been 
predicted. 

A set of three references (10,16,27) in Group 3 were used to predict 
the wind speed at sites (22-4) to determine if utilizing a reference group 
closer to the prediction group sites 22-24 would decrease the prediction 
errors. The delays for prediction of sites 22-24 utilizing references 
(10,16 ,27) were determined based on a 2 minute filtering interval rather 
than the 2 hour filtering interval used to determine the delays for 
reference group 1-5 given in Table 2. The delays for site 10, 16, and 27 to 
each to the prediction sites 22-24 are given in Table 5. The maximum delays 
to site 22, 23, 24 is 275, 280, and 290, respectively, using references 1-5 
but the average delay for sites 22-24 is 181, 253, and 221, respectively, 
using references 10, 16, 27. The delays decrease from 16-33% where the 
geographical distance between sites 22-24 to the reference sites has 
decreased approximately 40%. 
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Table 4. Table of wind speed prediction error for 2 minute, 10 min- 
ute, AND 2 HOUR DATA OF 05/02/79 (3 - 10 P.M.) 


SITE 

2 MINUTE 
FILTERED 
ERROR (MPH) 

10 MINUTE 
FILTERED 
ERROR (MPH) 

2 HOUR 
FILTERED 
ERROR (MPH) 

6 

5.20 

4.47 

1.51 

7 

5.32 

4.71 

0.51 

8 

2.91 

1.91 

0.49 

9 

5.52 

3.93 

0.23 

10 

4.65 

3.29 

1.00 

H 

5,82 

4.02 

0.28 

12 

4.83 

3.74 

0.48 

13 

4.28 

2.89 

0.31 

14 

4,78 

3.69 

0.64 

15 

3.49 

1.89 

0.11 

16 

4,64 

3r83 

1,06 

17 

3.53 

7.24 

0.66 

18 

7.98 

6.11 

0.36 

19 

4.81 

3.84 

0.75 

20 

2.80 

1.69 

0.14 

21 

4.02 

3.01 

0.47 

22 

8.12 

5.73 

0.36 

23 

3.97 

2.96 

0.15 

24 

8.59 

5.80 

0.33 

25 

4.74 

4.02 

0.68 

26 

5.02 

4.11 

1.22 

27 

4.24 

3.51 

0.67 
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Table 5 Comparison of delays using Sites 10, 16, 27 and 
Sites 1-5 as references. 



Table 6 Comparison of wind speed prediction errors using 
Sites 10, 16, 27 and Sites 1-5 as references. 
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The errors between the actual and pradictad wind speed records are 
given in Tfcbie 6 and decrease 0-25% using retaenoas (10,16,27) rather than 
(1-5) . The errors are for the individual site predictive nodal. 

Sites (7,9,11) are also used as references. These sites experience the 
s econ d storm but not the first. The delays given in Table 7 are comp u t e d 
with tbs 2 minute smoothed data and are considerably smaller than utilising 
references (10,16,27) or (1-5). These delays appear to rep res ent the delay 
asso c i a ted with the propagatio n speed on the wind shift (40 mph) and the 
distance (40 miles) between refere nc es (7,9,11) and prediction sites 22-24. 
The delay from reference 7 is longer then 9, 11, and is approximately 10-15 
miles further from the prediction sites. Site 24 is farther from the 
references (7,9,11) than sites 22 and 23 and again has longer delays. The 
use of a reference without the storm and a short 2 minute smoothing interval 
clearly and accurately determines delays associated with the wind shift 
propagation. The errors between the actual and predicted wind speed at 
sites 22-24 is 20-401 less than using references (7,9,11) rather than 
refe re nces (1-5) as shown in Table 8 for both the 2 and 10 minute filtered 
data. Figures 13 and 14 plot the actual and predicted wind speed record on 
sites 22-24 for reference (10,16,27) and (7,9,11) respectively. Note that 
the delays are much shorter using reference (7,9,11) and that the errors are 
very small except during the storm. A very surprising result is that the 
rise in wind speed and the cyclic variations in wind speed associated with 
the storm appear to be more accurately predicted on sites 22 and 23 using 
sites (7,9,11) that did not contain this first storm. The error during the 
storm is still quite large and the cyclic variations in the actual record 
are never truly predicted at any of the sites. 

These results further confirm the conclusion that one should not 
attempt to predict storms because 


(1) utilising references without the same storm and a short filtering 
interval captures wind shift propagation and produced substantially 
more accurate predicted wind speed records before, during, and after 
the storm; 

(2) the cyclic variation in the storm eithe. is quite different at 
different points in the storm or changes so radically over time that 
attempting to accurately predict the storm induced cyclic variation is 
impossible. 
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Prediction 
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1 

1 

j 

| 

Max torn 
Delay Fran 
Minutes 

^KTiSEh IHSIiSBil 

! 2 Hour Filtered 

2 Minute Filtered Data 

22 

1 

1 

i 

j 

275 

78 

33 

32 

23 

i 

! 

i 

295 

86 

43 

48 

24 

i 

i 

l 

290 

105 

67 

63. 


Table 7 Comparison of delay using Sites 7 , 9, 11 and 
Sites 1-5 as references. 



Reference 7, 9, 11 


Error toph) 

Error (roti) 

Error (npli) 


2 Minute Filtered 

10 Minute Filtered 


22 | 5.737 

3.123 

3.732 

5.463 

23 ' j 2.966 

3.973 

2.449 

3.230 

A | 5.802 

8.599 

4.885 

5.637 


Table 8 Comparison of wind speed prediction error using 
Sites 7 , 9, 11 and Sites 1-5 as references. 
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Figure 13. Comparison of actual and predicted wind speed for 3:00 - 10:00 p.m. 

USING 2 MINUTE MOVING AVERAGE FILTERED RECORDS AND SITES 10, 16, AND 
26 AS REFERENCES. 
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Figure 14. Comparison of actual and predicted wind speed for 3:00 - 10:00 p.m. 

USING 2 MINUTE MOVING AVERAGE FILTERED RECORDS AND SITES 7, 9, AND 
11 AS REFERENCE. 
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SECTION 5 


5.1 HIND ROWER SIMULATION AND PREDICTION 

Smoothing of wind speed records before simulation of power out of an 
individual wind turbine or wind turbine array is shown to cause serious 
distortion of the wind turbine or array power in Section 5.2. Smoothing of 
the wind speed records utilized to produce the predicted wind speed record 
and therefore cbe resultant simulated array power is shown to cause 
significantly less distortion from results presented in Section 53. The 
predicted wi nd speed and the predicted array power record that is produced 
from it are effected by an inherent spatial filtering caused by the use of 
several reference sites at different locations. Thus* predicted wind array 
power is not effected by the smoothing of the reference wind speed records 
since the spatial filtering has performed the smoothing. Distortion caused 
by the spatial filtering inherent in the wind speed prediction process is 
shown in Section 5.4 to cause significant wind array power prediction error# 
when the smoothing is performed on either the reference wind speed records 
used to produce the predicted wind speed or the actual wind speed record. The 
effects of utilizing two predicted wind speed records in the wind turbine 
cluster to simulate the wind power variation out of that wind turbine array is 
presented in Section 53. 

5.2 EFFECTS OF FILTERING ON THE SIMULATION OF WIND ROWER VARIATION 

The simulation of power out of a MOD-2 wind turbine given a record of 
wind speed utilizes a computer program developed in our earlier research 
[20]. The simulation program utilizes the nonlinear algebraic ptwer versus 
wind speed curve. The logic for high speed shutdown and startup and low 
speed shutdown and startup is included. The program can also simulate the 
power out of an array of MOD-2's given the siting configuration and the 
propagation speed of the meteorological event. The siting configuration is 
specified in terms of the abscissa and ordinate position in miles from the 
wind turbine that is at the boundary of the array in the direction of 
propagation. The square array simulated in this section is composed of 81 
sited in a matrix with separation of 1 mile in the abscissa and ordinate 
direction. The speed of propagation of the meteorological event is taken as 
30 mph# which is quite slow compared to the propagation speeds for 
meteorological events observed in our research on the SESAME array. A 
larger meteorological event propagation velocity increases the rate of 
change of power out of an individual wind turbine and out of an array of 
wind turbines. 

The effects of filtering the wind speed record on site 14 on the 
magnitude of the power out of a single wind turbine is shown in Figure 15. 

The wind speed is near rated velocity and thus the wind power is repeatedly in 
and out of saturation. The magnitude of the variation can be 13 megawatts 
and the periods of these oscillations range from 2 minutes to 10 minutes in 
the unfiltered case in Figure 15a. The 2 minute and 5 minute smoothing 
greatly reduce the magnitude of these variations so that the peak variation 
is 13 MW and 03 MW# respectively. The period of the oscillations increase 
also as the higher frequency components of these oscillations are filtered 
out. 


Figure ! 



















Figure 15c 
















The power out of individual wind turbines at site 14, 17, and 18 was 
simulated to indicate the variation in wind power for three wind turbines 
that are within 20 miles of each other as shown in Figure 16. The power out 
of the three wind turbines is given in Figure 17 utilizing wind records from 
1:00 - 6:00 pun. on May 2, 1979 with a 10 minute smoothing interval. The 
1.6 MW maximum wind power variation on site 14 in the unfiltered record is 
now only 0.5 MW and the period is now greater than 10 minutes. This 10 
minute smoothed power output record is not at all similar to that in Figure 
15a with no smoothing effects on the magnitude, frequency, or shape of the 
variation. Filtering the wind speed record before simulation of power 
totally destroys the power variations out of an individual wind turbine for 
a front. 

The wind speed increase and conoommitant wind power increase for the 
arrival of the front can be observed in sites 17 and 18 for the power 
simulated utilizing the 10 minute smoothed records. The magnitude of the 
power output is very different because the average and rate of change of 
wind speed is quite different on these three sites. Site 17 eventually 
nears and exceeds rated wind speed on a MOD-2 at t ■ 230 minutes but site 18 
wind speed and power never reach rated velocity at t - 300 minutes. The 
shape of the wind power variation at sites 17 and 18 are very different even 
though the majority of the variation has been filtered out with the 10 
minute smoothing interval. 

Thble 9 shows the mean and rms errors between the filtered and 
unfiltered individual wind turbine power output. Note as the smoothing 
interval increases, from 5 to 30 the mean and rms error in power changes 
relatively little for site 14 but shows larger changes for site 17. The m + 
3 error is .85 MW for site 14 and 1.35 MW for site 17 with the 5 minute 
smoothing interval. This error is very large compared to the 2.5 MW 
capacity of the fOD-2 wind turbine. 

The power output of the 81 wind turbine array at site 14, 17, and 18 is 
given in Figure 18. The 10 minute smoothed individual wind turbine power 
record is multiplied by the number of wind turbines in an echelon and 
delayed by 


d^ - distance of echelon i from the first wind turbine in the array in the 
direction that the meteorological event is propagating 

V 0 - speed of propagation of the meteorological event 

to produce the power out of the ith echelon of wind turbines. The power out 
of the array is the sum of the power out of all echelons. The power out of 
the 81 wind turbine array is sites 14, 17, and 18 shown in Figure 18 is a 
smoothed version of that in Figure 16 for the individual wind turbine. The 
array simulation procedure is a spa tied filtering process and thus the 
maximum percentage change in power variation is only 10% on the site 14 wind 
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Figure ; 
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Table 9. Table of error for actual wind power on a single wind 
TURBINE FOR FILTERING THE WIND SPEED WITH 5 , 10, AND 30 
MINUTE SMOOTHING INTERVALS. 


HMI 

f'EAN ERROR BETWEEN 5 MINUTE 
FILTERED AND UNFILTEREE DATA 

RMS ERROR BEIVEEH 5 MINUTE 
.FILTERED AND UNFILTERED DATA 

.14 

.13 

.24 

17 

.27 

.36 

SITE NO. 

MEAN ERROR BETWEEN 10 MINUTE 
FILTERED AID UNFILTERED DATA 

RMS ERROR BETWEEN 10 MINUTE 
FILTERED AND UNFILTERED DATA 

14 

0.132 

0.247 

17 

0.233 

0.371 


WfaMaiEH 


14 

0.141 

0.272 

17 

0.317 

0.419 
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Figure 18 . Actual wind power from an 81 wind turbine array sited near Sites 14 

17 J AND 18 USING 10 MINUTE MOVING AVERAGE FILTERED DATA. 



turbine array but 30% from tba ait* 14 wind turbine. The period of tba 
variations has incraaaad for tba arrsy power output dua to tbia apatial 
fil taring. Sita 17 approaebaa but navar raacfaaa ratad array powar output 
although tba individual wind turbina raaebad ratad powar output at t • 210 
minutes. Similarly, aita 18 navar avan approaebaa ratad array powar output 
although tba individual wind turbina powar approachad ratad powar at t ■ 
300. 


Oonparlson of tba array powar utilizing an unfiltered, and a filtered 
wind speed record ia shown in Figure 19a-d for 2, 5, 10, and 30 ainuta 
smoot*»trjg intervals. Tha filter with 2 ainuta smoothing interval appear to 
eliminate tba cyclic turbulancr*variation and reduoa total powar output by 
4%. The filtering with 5, 10, and 30 ainuta smoothing intervals 
significantly reduced tba powar produced from tba array although tba ahape 
of the array variation ia retained until a 30 ainuta smoothing interval ia 
used. Tha percentage reduction in powar out of tba array can reach 20% due 
to this filtering but raaebad greater than 70% for a single wind turbine. 

The moan and rms error between array powar output baaed on the 
unfiltered and a filtered wind speed record ia given in Table 10. Tba mean 
and rms values increase with smoothing interval and are less than 81 times 
the values for an individual wind turbine in Table 9. This error increases 
with smoothing interval more than for the individual wind turbine. The m + 
3o magnitude of 32.1 for site 14 and 41 for site 17 and a five minute 
interval which ia still a significant percentage of the 203 MW array 
capacity. 

5.3 EFFECTS OF FILTERING ON PREDICTED ARRAY SOWER 

The previous section investigated the effect of filtering of an actual 
wind speed record (at a meteorological tower site in the SESAME array) on the 
power out of a fictitious 81 turbine array sited on a 9 mile x 9 mile square 
area behind the meteorological tower. It was found that filtering 
significantly increased the power of an individual wind turbine and out of 
an array and significantly reduced the cyclic variation out of a turbine and 
out of an array. 

Prediction of wind speed is a spatial filtering process where the wind 
speed at several sites are weighted and delayed and then summed to produce 
the predicted wind speed record. Thus, filtering these wind speed records 
before the predicted wind speed record and the predicted power out of the 
array is com p u t ed will be shown to have very little effect compared to the 
filtering of the actual wind speed record which is than used to simulate the 
power out of an array. The reduction of the effects of filtering the 
reference wind speed records on the predicted wind power variation out of 
the array is due to inherent spatial filtering in the wind sp ee d prediction 
process. This spatial filtering of the wind speed prediction will be shown 
to incresse the power predicted out of the array compared to that produced 
from the actual wind speed record. This optimism in predicting the trend 
could reduce the spinning and operating reserve provided by the unit 
commitment and thus the reliability of the power system if wind prediction 

error Q^(k ♦ j/4 - D were set to zero in the unit commitment procedure 
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TIME IN MINUTES 


Table 10 . Table of error for actual array power for filtering 
THE WIND SPEED WITH 5j 10 , AND 30 MINUTE SMOOTHING 
INTERVALS. 



SITE NO. 

MEAN ERROR BETWEEN 5 MINUTE • 
FILTERED AND UNFILTERED DATA 

RMS ERROR BETWEEN 5 MINUTE 
FILTERED AND UNFILTERED DATA 

«* 

6.10 

3.72 

17 

8.11 

10.98 

SITE NO. 

MEAN ERROR BERCEN 10 MINUTE 
FILTERED AND UNFILTERED DATA 

RMS ERROR BETWEEN 10 MINUTE 
FILTERED AND UNFILTERED DATA 

14 

7.768 

11.791 

17 

9.037 • 

11.124 

SITE NO. 

MEAN ERROR BETVEEN 30 MINUTE 
FILTERED AND Uf FILTERED DATA 

RMS ERROR BETWEEN 30 MINUTE 
FILTERED AND UNFILTERED DATA 

14 

9.307 

14.807 

17 

13.471 

16.153 
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developed in Section 2. Q»(k + j/4 - 1) should not be set to zero but 
should be based on an estimte on the maximum error in the trend wind power 
prediction Thus, this optimism in predicting trend wind power should not 
affect operating reliability if QjJ(k+j/4-l) is set proportional to the wind 
power prediction error. 

Ihe inherent spatial filtering in the wind speed prediction process 
will also be shown to eliminate the cyclic variation in the actual wind 
power variation. 

The predicted power out of an array at site 14 utilizing filtered and 
unfiltered reference wind speed records for wind speed prediction is shown 
in Figure 20a-c for 2, 5, and 10 minute smoothing intervals. The difference 
between the array power utilizing filtered and unfiltered records at these 
reference sites for wind speed prediction and thus predicted wind array 
power variation is negligible as shown in Figures 20a-c. 

The mean and rms error between unfiltered and filtered wind array power 
variations are tabulated in Ihhle 11 for actual and predicted wind power 
variation Ihe mean and rms error for predicted power at sites 14, 17, 18 
are less than for the actual wind power due to inherent filtering provided 
by the spatial filtering in the prediction process. 

5.4 COMPARISON OF PREDICTED AND ACTUAL WIND ARRAY POWER 

Ihe previous two sections have shown that the effect of filtering the 
wind speed record at a site used to simulate the actual wind array power is 
large. Filtering the reference wind speed record used to predict the wind 
speed at that site, which is then used to simulate predicted array power, was 
shown to have little effect on predicted array power variation. A 
comparison of actual and predicted wind power variation must thus be made 
when no filtering is performed on the actual wind speed record at the site 
or on the reference wind speed records used to predict wind speed at the 
site. 


A comparison of the actual and predicted wind array power at site 14 is 
plotted in Figure 21. The mean and rms error is 5.78 MW and 9.70 MW 
respectively. This plot shows the high frequency cyclic variation in the 
actual which is not present in the predicted wind power variation due to the 
spatial filtering of the wind speed prediction process. The predicted wind 
power record is consistently greater than the actual wind power record as 
expected due to the inherent spatial filtering in the wind speed predictor, 
ihe m + 3a error * 34.5 MW which is approximately 15% of the array capacity. 
Use of only one measurement site to simulate wind array power variation can 
cause additional error as shown in the next subsection and in Section 6 of 
this report. 


Figure 20a. Comparison of predicted wind array power at Site lf» using 2 minute 

MOVING AVERAGE FILTERED AND UNFILTERED DATA. 
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Figure 20b. Comparison of predicted wind array power at Site M using 5 minute 

MOVING AVERAGE FILTERED AND UNFILTERED DATA. 




Figure 20c. Comparison of predicted wind array power at Site M using 10 minute 

MOVING AVERAGE FILTERED AND UNFILTERED DATA. 
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Table 11 . Error in actual and predicted array power variation for 

UNFILTERED AND 2 t 5j AND 10 MINUTE FILTERED WIND MEASURE- 
MENT RECORDS. 
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5.5 EFFECT OF SEVERAL WIND SflSBD PREDICTION RECORDS ON WIND ARRAY POWER 

PREDICTIONS 

The wind turbine array are generally spread out over a large 
geographical area due to the fact that one is cautioned to site wind 
turbines no closer than ten blade diameters apart. Although wind turbines 
may be sited more closely than this guideline? not all sites may have 
excellent wind resource and not all of the sites with excellent resources 
may be available due to other uses? and environmental concerns. Thus? a one 
mile on a side matrix sited away may be typical for seme areas. 

The 81 wind turbine array with a matrix one mile on a side separation 
between wind turbine covers a 9 mile by 9 mile area. Our research on wind 
speed prediction in Section 4 has shown that the average and rate of change 
of wind speed at sites 10 miles apart can be quite different. Thus? it 
appears that one wind measurement site will not be sufficient to adequately 
predict wind speed and power at all wind turbines in an array covering 81 
square miles due to the fact that average wind speed? wind speed variation? 
and rate of change of wind speed varies widely over the array and no attempt 
is made to assess this geographical variation if one wind measurement site 
is used. 

The research performed in this section is to confirm the need for more 
than one wind speed prediction site in each wind turbine cluster. Wind 
power will be predicted for an 81 wind turbine array sited at Sites 14 and 
17 as shown in Figure 22. The predicted wind speed at Sites 12 and 14 will 
be utilized to simulate array power at the array at Site 14. The predicted 
wind speeds at Sites 15 and 18 will be utilized to simulate the power out of 
the array at Site 18. Note the measurements at each array are at the front 
and back of the array in the direction of meteorological event propagation 
for the front observed from 1:00 - 6:00 ptm. on Hay 2? 1979 on these sites. 
Sites 23 and 25 are utilized aw reference measurement sites to predict wind 
speed at 15? 18? 12? and 14 rather than reference Sites 19? 22? 23? and 25 
utilized in the previous work. The wind speed record utilized to predict 
wind speed and thus simulate predicted power or the actual wind speed used 
to simulate actual power are filtered with a 2 minute smoothing interval for 
all cases in this section. 

The actual and predicted array power at Site 14 utilizing only the 
actual and predicted wind speed at Site 14 is shown in Figure 23. The mean 
and rms errors are 6.93 and 10.93 but the predicted wind power is less than 
the actual wind power record. This result may be explained by the fact that 
only two reference sites are utilized that are very close together and have 
similar wind speed records. The spatial filtering of the wind speed 
prediction process is thus quite limited and explains why the predicted 
power is not larger than the actual power. The m + 3 a error of 40 is 20% 
of the total array capacity. 

The actual power out of the array at Site 14 utilizing solely the 
measurement record at Site 14 along with the predicted array power utilizing 
the average of the predicted wind speed at Site 12 and 14 is plotted in 
Figure 24. The mean and rms error is 5.99 and 8.76 but the predicted power 
is now larger than the actual array power. The predicted wind speed at Site 
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Figure 22. Map of wind turbine array siting for the wind turbine array that 
utilizes measurements at Sites 12 and 14 and for the array that 
utilizes measurements at Sites 15 and 18. 


Figure 23 . Comparison of actual and predicted wind power that utilizes wind records 
at Site 14 . 
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12 is larger than at Sib a 14 raising the power out of the array by an 
average of 12 megawatts. The shape of the predicted array power utilizing 
two wind prediction sites is closer to the actual power based on only Site 
14. h comparison of the actual and predicted power utilizing both an 
average of the actual and predicted wind speeds was not carried out on this 
data but was carried out on the Site 18 array. 

The power out of the array using the actual and predicted wind speed at 
Site 18 is shown in Figure 25. The predicted wind power is greater than the 
actual wind power as expected due to the spatial filterirg associated with the 
prediction process. The mean and ms error is 19J.4 and 22.67 which is very 
large. The error m ♦ 3 o * 87 MW which is 40% of the array capacity. 

Although these errors are larger Figure 24 indicates the predicted power 
tracks the actual power record ramp increase but the site specific cyclic 
variations with a 30 minute period are not predicted. The actual and pre- 
dicted power from the array based on an average of the actual and predicted 
wind speed records at Sites 15 and 18 is plotted in Figure 26. Die pre- 
dicted power is again larger than the actual wind power from the array. The 
mean and rms errors are 18.41 and 21.94 and the m ♦ 3 a error - 84.5 is 
large. The most important difference utilizing the two Sites 15 and 18 
rather than just Site 18 is very substantial change in power out of the 
array. The power out of the array may be as much as 35 MW more utilizing 
the average wind speed at 2 sites than utilizing one site. A substantial 
increase in power output utilizing the average wind speed at 2 sites com- 
pared to utilizing just a single measurement or predicted wind speed was 
also observed on the Site 14 array. The wind speed at Site 15 was larger 
and had different site specific variations than at Site 18 and thus the use 
of two sites reduced the magnitude of these site specific effects of both 
the actual and predicted wind array power records utilizing the average of 
the wind speed at two site:}. 
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Figure 25. Comparison of actual and predicted wind array power using wind records 
at Site 18. 
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Figure 26. Comparison of actual and predicted array power based on an average 
OF WIND RECORDS AT SITES 18 AND 15. 
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SECTION 6 


COMPARISON AND ACCURACY OF 
HIND ARRAY POWER SIMULATION (CTBODS 


lbs results of the previous section indicate that utilizing several 
reference sites to predict wind speeds at one or more sites in a wind turbine 
array reduces site specific effects of any specific reference site but 
effectively performs a spatial filtering of the wind speed profile associated 
with the specific meteorological event. The spatial filtering reduces the 
level of the power predicted from an array in a manner similar to time 
filtering of the predicted wind speed records. This negative bias in the 
predicted wind power is unavoidable if several reference sites are utilized to 
predict wind speed at sites in the array. 

A preliminary study of the effects of utilizing several prediction sites 
to simulate the wind power from an array of wind turbines was conducted. The 
results suggested that each prediction site can have very different mean r 
standard deviation/ and site specific cyclic variations even though the sites 
may be within five miles of each other. These site specific effects 
greatly effect wind array power if only one site is used but the site 
specific effects are reduced if several sites are used. It is clear that 
utilizing as many prediction sites as possible would increase array power 
prediction accuracy. Since the number of prediction sites is limited/ due 
to cost of instrumentation and communication hardware for measurements at 
each additional prediction site and the computational requirements for 
producing each additional wind speed prediction record/ care must be uti- 
lized in simulating array power variation to minimize the site specific 
effects and minimize the wind array power prediction error from a limited 
nunber of wind prediction sites. 


The results thus indicate a more complete study is necessary to assess 
the accuracy and the differences between the array power obtained using six 
different methods of simulating wind array power variation. These different 
simulation methods utilize one or several wind prediction sites to simulate 
array power variation. 


Figure 27 shows a wind turbine array composed of N wind turbines each 
sited at a distance {d^ N from the reference wind turbine in the array that is 

0*1 

the first wind turbine to be effected by the wind variation of the 
meteorological event. The J wind measurement sites are sited at a distance 
from the reference wind turbine in the array. The wind speed record at site 
*Mt) is reflected to occur at the reference wind turbine site by advancing 
the record by T j ■ Dj/ to produce Wj(t + Tj) where V Q is the speed of 

V 0 

propagation of the meteorological event. 


'3 


The predicted wind power from the set of N wind turbine based on the 
predicted wind speed record is 
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Figure 27. Wind turbine array configuration showing distances Dj 

FROM WIND PREDICTION SITES AND DISTANCE D*. TO WIND TURBINE 



P M (t) - “ Fj (t + Tj - y 
n-i 

where the power from the reference wind turbine is 
Pj (t + Tj) - f{Wj(t + Tj) } 

and where f( * > is the nonlinear algebraic relationship between wind power 
versus wind speed for a particular wind turbine model. The wind array power 
array (t) is the sum of the power from each of the individual wind turbines 
in the array Pj(t + L * L) produced by delaying the simulated wind power 

J J ^ 

P-:(t + Tj) at the reference wind turbine site where = -- . 

J J v o 

The second method averages the wind speed record reflected to the site of 

the reference wind turbine record 


w A (t) 


1 J 

% E Wj(t + Tj) 
J j*l J J 


The power produced by a wind turbine from this average wind speed record is 
P A (t) = f{W A (t)} 

She power from the array of wind is then produced by delaying and summing the 
averaged reference wind turbine power record P A (t) for each wind turbine in 
the array 

p AII (tJ = 2 p A (t ~ V 
n=l 


N i J 

= Z f(r Z Wj(t + Tj - t-)) 
n=l J j=l J J 

The third method simulates the power from each of the wind speed records 
reflected to the site of the first wind turbine to be effected by the 
meteorological event. 

Pj (t + Tj) = f(Wj(t + Tj)) 

An average wind power record 

1 J 

= - z Pj(t + Tj) 

J j*l J J 


e A (t) 
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is then produced. Hie array power is obtained by delaying and summing the 
average reference wind turbine record for each wind turbine in the array 


N * 

P AIII (t) * Z . P A tfc “ ^ 
n»l 

1 N J 

* r Z Z f(W(t + Tj - U) 
J n-1 j-1 2 


The fouttb method determines the power from the subset of wind turbines 
closest to each wind prediction site j based on the wind speed prediction 
record JMt). The wind measurement sites j are assumed to be ordered based on 
distance "bj > Dj_^ from the reference wind turbine site. Similarly f the 
wind turbine sites n are assumed to be ordered based on the distances d- > 
from the reference wind turbine site. The jth subarray of wind turbines 
utilize wind prediction record Wj (t) are located in the interval 


D 3-1 +D j 


Id i 


D f D j+l 

2 


as shown in Figure 27. Defining 


j=l 


N(j) 


nexCn/c^ < 


°i +D j + i 
‘ ~2 


} 


j / 1 


the wind turbines n belonging to the jth subarray satisfy 
N( j) < n 1 N(j + 1) 


The power from the jth subarray is 
I N 



p sj it) 


where 


< 


N( j) 

Z Pj(t + Tj - 
n=N(j-iy+l 2 




N(l) 

Z 

n=l 


*i(t + T i “ ^ 


j s J 

j = 2 f 3/..» (iW 
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f {Wj (t) } 

. htsZh 




1*1 


j - 1 


( 12 ) 


T 3 


D i +D i-1 

o r a r i~) 



3*1 


( 13 ) 


j ■ 1 


The values of foe each array are chosen so that wind measurement site j 

D j +D j-1 

within subarray j is reflected to the boundary ----- — between subarrays j-1 

and j. The wind turbine power record Pj(t + T-) at the boundary between 
subarray j-1 and j is then delayed by t^ based J on the distance of the wind 
turbine site n in subarray j from the boundary between subarrays j-1 and j and 
summed for each wind turbine site in that subarray to produce P s j(t). The 
power out of the J subarrays that constitute the array is the sunr of the 
subarray power 


J 

p Alv(t) ■ ^ p sj (t) 

The fifth method of simulating power out of the array again simulates 
power out of each subarray separately and then sums the power out of all 
subarrays to produce the power out of the array. The jth subarray of wind 
turbines utilize an averaged wind prediction record 


and 


W. 


jj+1 


(t) = - [IMt) + W,(t + T. 


3+1 


cover the same subarea defined for 


'j-l- u j . , . J j+1 +U j 
"~2 1 d 1 2 


- Tj) ] 

the previous method 


The wind speed prediction record w -H+i(t) is the average of the record Wj(t) 
and the record Wj + ^(t) advanced by Tj^ - Tj to occur at site j. J 

The power produced based on this averaged wind speed record for subarray 

j is 


P jj+1 (t) - f(W jj+1 (t» 


The power from the jth subarray is obtained using the same method used in 
method 4 where Pj(t) is replaced by Pjj+i(t) and thus 


N(l) 

I P! 2 (t + T x - t^) 


n«l 


p sjtt) 


N( j) 

Z 

n*N(j-l)+l 


p jj+l (t + T j - ^ 


j*2 /3 1 • • » 


N(J) 


Z 

np«N(J-l)+l 


p J-l,J (t + T J “ ^ 


where t n and T^are defined in Q.2) and CL3). The power out of the array is 
again the sum of the power out of the subarrays 



J 

z P gi (t) 
j-i 3 


The sixth method of simulating the power out of an array is very similar 
to the fifth method. The average power for sites in the jth subarray is 


- 2 - i p j 


where 

Pj(t) = f(Wj(t)) 

Pj+l(t + Tj+i - Tj) * f(Wj(t + Tj + ^ - Tj)) 

The power based on the wind prediction site j+1 is advanced by T i+1 - to 
appear as if it occurs at site j. J J 
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Figure 28. Vertical wind turbine array configuration showing 
LOCATION OF WIND PREDICTION AT SlTE 22, 18, AND 
15 AND THE 90 WIND TURBINES IN THE ARRAY. 
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The power from subarray j Is then 


N(j) 

t + T i ' tn) 

N(l) 

E «i2(t + *1 - y 

n-1 

N 

I P* T I 7(t + T-t ** t-.) 

n-H(J-l)+l ' W,J J ^ 

The power out of the array is the sum of the power out of the subarrays 


p sj 


p array^ * ^ p sj ^ 

The power from an array of 90 wind turbines sited in a rectangular area 
of 9 miles by 9 miles r is simulated using the six different methods just 
described. The spacing between wind turbines in both the latitudinal and 
longitudinal direction is one mile. The wind prediction sites are located 
within this rectangular an: ray at sites 22, 18, and 15 as shown in Figure 28. 
The prediction of the wind speed at sites 22, 18, and 15 is based on wind 
speed measurements at sites 19 and 23. 

The wind power simulated based on the single wind speed measurement at 
site 22 is shown in Figure 29 along with the power produced based on methods 2 
and 3. Methods 2 and 3 produce an average wind turbine based on either 
averaging the wind speed or wind turbine power records from sites 18, 22, and 
15 reflected to the reference wind turbine site. The array power is then 
produced in both methods by delaying and summing this average wind turbine 
power record for each wind turbine site in the array. These average wind 
array power methods 2 and 3 using wind prediction sites 18, 22, and 15 produce 
wind turbine array power records that are smaller in total power than the wind 
power record based on method 1 using the single wind speed record at site 22. 
This result occurred because the magnitude of the wind speed measured at sites 
18 and 15 are smaller than the wind speed record measured at site 22. The 
array power produced based on the single wind speed measurement at 18 and 15 
is shown in Figure 29b and 29c, respectively, along with the averaged wind 
power record based method 2 and 3 using prediction sites 15, 18, and 22. 

These results show that the wind array power simulated based on single wind 
speed measurement at either 15 or 18 is approximately 30% less than that 
produced based on the average of 15, 18, and 22 and approximately half of that 
produced f r om the single wind speed measurement at 22. Since all three sites 
are located in the 90 square mile geographical area containing the 90 wind 
turbine sites in the array, it is clear that wind turbine array power 
simulation based on a single wind measurement record correctly captures the 
shape of the trend change seen in the array but does not accurately capture 
the magnitude of the power produced by averaging the wind speed or wind 
turbine power records before simulating wind array power. The very large 
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cyclic variations that are measurement site dependent are drastically reduced 
fay averaging the wind speed records as indicated in Figure 29 a, br and c. 

The wind array power based on method 3 is virtually identical to that 
produced by method 2. This result suggests that the array power produced by 
first averaging wind speed before simulating the power from an "average wind 
turbine" gives nearly identical results to simulating power for each wind 
speed record and then averaging to produce the "average wind turbine 
record." The array power is produced by delaying the average wind turbine 
record for each wind turbine and then summing to produce the wind array 
power record in both methods. 

It is clear from the results shown in Figure 29 a, br c that utilizing 
several wind measurement widely dispersed geographically within the wind 
turbine array is necessary to accurately simulate wind array power. The 
wind array power simulation results for methods 4, 5, and 6 should indicate 
whether 

(1) simulating power from subarrays based solely on the wind measurement site 
or sites closest to a subarray of wind turbines and then summing to 
pcoduoe a total array wind power record will result in significant 
differences from averaging wind speed (method 2) or power (method 3) and 
then simulating the total array wind power; 

(2) whether use of closest single measurement site or the average of the 
closest two wind measurement sites can effect the power produced from a 
wind array power simulation method. 

The results to be presented will show both of the above factors can cause 
significant changes in the power simulated from an array of wind turbines when 
the wind speed measurement records at various sites are quite different. The 
power produced from the array using method 4, 5, and 6 are shown in Figure 
30a, b, and c respectively, along with the power simulated from the array 
using methods 2 and 3. The power produced from the array using method 4 is 
considerably smaller than for methods 5 and 6 even though the number of wind 
turbines in the subarrays 1, 2, and 3 for the three methods are identical. 

Wind speed record 22 is utilized solely in subarray 1 in all three methods and 
thus the power out of the subarray is identical in all three methods as shown 
in Figure 31a. An average wind speed record (14) based on W^g(t) and ^22^ i® 
utilized to produce subarray power in method 5 and an average of the powers 
produced based Wi g(t) and W 22 (t) is used to produce the subarray power for the 
second subarray in method 6. The subarray power for the second subarray 
depends solely on W^g(t) for method 4. Since W 22 <t) is much larger than 
Wig(t), the power produced by method 5 and 6 for subarray 2 are much larger 
than for method 4 as shown in Figure 31b. 

The power produced for subarray 3 is very similar for all three 
simulation methods as shown in Figure 3.1c. Method 4 utilizes W^c(t) alone to 
simulate power from the array. Method 5 averages Wig(t) and W^ft) to produce 
power In the third subarray and method 6 averages the power produced based on 
W^e(t) and W^g(t). It is clear that since the wind speed records are similar 
methods 4, 5, and 6 give nearly identical results in subarray 3 from the 
results in Figure 31c. 
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The effect of siting wind turbines in different patterns was investi- 
gated. The 90 wind turbines were sited in a horizontal pattern as shown in 
Figure 32 in addition to the vertical pattern shown in Figure 28. The 
differences in the array power produced by methods 2-6 for the vertical and 
horizontal array configurations is shown in Figures 29 and 33 respectively. 
Mote as the number of wind carbines in a straight line normal to propagation 
of the front increase* the size of the cyclic variations increase the ramp 
rate of change of power appears to be quite similar on all three siting 
configurations. The differences between arrays over the simulation methods 2- 
6 on the vertical array configuration are the same as the differences observed 
between the simulation methods 2-6 for the horizontal array configuration. 

One should attempt to minimize the impact of any one single wind 
measurement on the power simulated from the array since one cannot be sure of 
how many wind turbines in an array will have similar wind speeds as a 
particular wind speed measurement site and one cannot be certain the wind 
speed measurements at any site are not site specific. In order to minimize 
the error between the power produced from simulation and power actually 
produced from the array* each measurement sites effect on total array power 
should be equal. Thus* the number of wind turbines in each subarray should be 
equal. Either wind speed at only one site or the wind speed measurement at 
the closest two sites may be averaged to simulate power in a subarray if the 
number of wind turbines in the subarrays are equal. 
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Figure 32. Horizontal turbine array configuration showing 

LOCATION OF WIND PREDICTION SITES AND THE 90 
WIND TURBINE ARRAY SITES. 


150 


Figure 33 . Comparison of simulated wind power variation of Methods 4 , 5 , and 6 

WITH POWER SIMULATED USING METHODS 1 , 2 , AND 3 . 




SECTION 7 


WIND POWER ERROR PREDICTION (GTBODS 

The wind power prediction error will generally exceed 15% of the wind 
generation level out of the wind turbine array. The wind power prediction 
error can be as large as the total array generation capacity if low speed and 
high speed shutdowns of the wind turbines in the array are not predicted. 

The large storm induced cyclic variations cannot be predicted using the 
present methodology and have been shown to cause repeated low speed shutdown 
and startups that very quickly reach wind turbine generation capacity. High 
speed shutdown may not always be predicted since the wind power prediction 
will always overestimate the actual wind array power level due to the 
spatial filtering caused by utilizing several reference wind measurement 
sites. 


The error in predicting wind power will normally be approximately 20% 
even when no storm related events occur. This wind power prediction error is 
associated with error in predicting wind speed at severed sites in the array 
and the error due to use of only a few wind prediction sites to simulate power 
from the large number of wind turbines in the array. If the wind power 
prediction error was 2% or less as it is for predicting load power variation 
in a utility/ this wind power prediction error could be neglected. The wind 
power under open loop control/ where the wind turbine array power output is 
not controlled/ is considered as negative load. The wind array power produced 
under open loop control reduces the load power in a utility that must be 
served by non-wind generating units and is not controlled just as load is not 
controlled and thus acts as negative load. The large error in predicting wind 
power would thus increases the effective error in predicting the load power 
that must be served by conventional generation. Operating reserve/ spinning 
reserve/ unloadable generation reserve/ and load following reserve levels 
within unit commitment schedules/ that determine when generating units must be 
brought on-line or shutdown/ must be adjusted based on the magnitude of the 
wind power prediction error. Modified generation control methods must be 
employed to compensate for the large cyclic wind power variations that 
comprise a portion of this wind power prediction error. Modified unit 
commitment and generation control methods that utilize this wind power 
prediction error are discussed in Section 2 of this report. 

Methods for estimating or predicting the error in wind power prediction 
are discussed in this section. A band on the wind power prediction is 
required since sufficient spinning reserve on conventional generation is 
required to cover cases where the predicted wind generation is greater than 
actual wind power generation and sufficient unloadable generation reserve is 
required when predicted wind generation is less than actual wind power 
generation. 

The first method to estimate wind speed prediction error is based on the 
assumption that the wind speed prediction error is time invariant^and is 
normally distributed. If the wind speed prediction error W(t) - W(t) is time 
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invariant and normally distributed, with zero mean and standard deviation w , 
then 

P{ W(t) - W(t) < 3c w } ~ .99 

A method for estimating wind power prediction error given the wind speed 
is normally distributed and the correlation between wind speeds records at 
different wind turbine sites is given in [2] and is discussed in Section 3 of 
this report. This method of predicting when power variation could be 
applied based on the following tests showing wind speed prediction error is 
time invariant and normally disturbed. 

The actual wind speed and predicted wind speed records for sites 21, 22, 
24, and 26 during the movement of a front from 1-6 pm. on May 2, 1979 is 
given in Figure 34 and 35. The actual wind speed records 23 and 25, used to 
predict wind speeds records at sites 21, 22 , 24, and 26 are filtered using a 
10 minute moving average filter in Figure 34 and by a 2 minute moving 
average filter in Figure 35. The rms errors for all sites is given in Thble 
12 for both the 2 minute and 10 minute filtered records. The errors 
decrease as the filtering interval increases. This result is confirmed by 
noting the error magnitude between the actual and predicted wind speed 
records for the 10 minute and 2 minute filtered records in Figures 34 and 35, 
respectively. The actual wind speed prediction error is plotted in Figure 
36 and shows the error has zero mean and appears time invariant. Thus, 
these two assumptions required for estimating a band around the predicted 
wind speed W(t), which will contain W(t) 

W(t) + 3c w iW(t) < W(t) + 3c w 

are satisfied. A statistical hypothesis test, that determines whether the 
wind speed error sampled every minute is normally distribution, was applied. 
The hypothesis test was applied to each SESAME array where wind speed 
prediction was attempted. The hypothesis that the wind speed prediction 
error is normal was accepted using a threshold on differences between the 
distribution of the sampled error and actual normal distribution. This 
threshold was selected so that 95% of such decisions would be correct and 5% 
of such decisions would be incorrect. The results of the statistical tests 
on every wind speed prediction site is given in Table 13 for both the 10 
minute filtered and the 2 minute filtered records. 

The wind speed prediction errors was decided to be normal in 16 of the 25 
sites when the 2 minute filtered data was tested for normality. Ten of 25 
sites were decided to be normal when the 10 minute moving average filtered 
data was used. The filtering eliminates and distorts the error variations 
causing more sites to reject the normality hypothesis. Sites 27, 17, 16, 21, 
19, and 6, where the normality hypothesis was rejected, are located at the 
eastern and western of the SESAME array where wind speed prediction using 
sites 23 and 32 was not as effective. Sites 8, and 12 had site specific 
effects that could have caused the error to be non-normal, which explains why 
the normality hypothesis was rejected for these sites. 

The results on this May 2nd 1-6 pm. data were very encouraging since the 
wind speed prediction error was shown to be normal, zero mean, and time 
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Table 12. Table of prediction errors for 2 minute and 10 minute 
FILTERED (1“6 P.M.) DATA OF MAY 2 , 1979 WITH SITES 
23 AND 25 AS REFERENCE. 


SITE no. ERROR (ilPH) ERROR (M°H) 

2 MIMUTE FILTERED 10 MINUTE FILTERED 



Figure 36. Error in wind speed prediction at Sites 21, 22, 24, and 26. 
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Table 13. Results of hypothesis test for prediction error at 2 and 
10 MINUTE FILTERED 1-6 P.M., DATA. SET OF MAY l , 1979. 


SITS NO. 

10 MINUTE FILTERED 

2 MINUTE FILTERED 

1 

reject 

reject 

2 

reject 

accept 

— * 
j 

reject 

accept 

4 

Gcceot 

accept 

5 

reject 

accept 

6 

reject 

reject 

7 

reject 

accept 

3 

reject 

reject 

9 

reject 

accept 

10 

reject 

accept 

11 

accept 

accept 

12 

reject 

reject 

13 

accept 

accept 

14 

reject 

reject 

15 

accept 

accept 

16 

reject 

reject 

17 

reject 

reject 

13 

accept 

accept 

19 

accept 

reject 

20 

accept 

accept 

21 

reject 

reject 

22 

accept 

accept 

24 

accept 

accept 

26 

accept 

accept 

27 

reject 

reject 
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invariant when the wind speed prediction is being successfully accomplished. 
The error band on the actual wind speed is plotted along with the actual wind 
speed in Figure 37. It clearly indicates the actual wind remains in the band 

W(t) » 3c w i W(t) 1 W(t) + 3a w 

and that the band is not too small or large based on the fluctuations of W(t) 
within the band. 

This method of estimating wind speed prediction error was retested for 
the wind speed prediction for wind records taken from 3 to 10 pun. on Hay 2nd. 
These wind speed records contained severe cyclic storm induced wind 
variation. Reference sites 7, 9, and 11 were chosen as references since they 
did not contain the storm induced variations, thing reference sites without 
storm induced cyclic wind variation was found to much more accurately predict 
the trend change in wind speed associated with a storm front. Since the 
cyclic storm induced variations are not correlated between sites due to their 
time varying characteristics, no effort is made to predict these variations 
when reference sites are chosen so that they do not contain such variations. 
The wind speed prediction error estimation procedure should hopefully be 
capable of estimating the magnitude of these variations. 

The error between the actual and predicted wind speeds for the sites 
where wind speed was predicted on the 3 to 10 p.m. May 2, 1979 wind record 
using references 7, 9, and 11 are given in liable 14. The errors are 
considerably larger due to the large cyclic variations due to the storm that 
are not predicted. The error magnitude was once again smaller if the 10 
minute moving average filter rather than a 2 minute, moving average filter is 
used to smooth the actual wind speed at each site and the reference wind speed 
records used to produce the predicted wind speed records. The hypothesis 
test was again applied to determine if the error between the actual and 
predicted wind speed records at each site was normally distributed. The 
hypothesis given in Table 15 that the error was nortnally distributed was 
accepted on only 7 of 24 sites for both the 2 minute and 10 minute moving 
average filtered records The sites where the hypothesis was accepted were 
generally close to the reference sites. The large cyclic variation due to 
the storm can be seen in Figure 38 where the actual wind speed record W(t) 
and the upper W(t) + 3 a w and lower W(t) - 3o w limit are plotted for sites 
22-25. Note that the cyclic variation that occurs between 300 and 380 
minutes are quite different in shape and magnitude at these different sites. 
Thus, the period of the cyclic variation due to the storm was cut out of the 
wind speed record for each site and the hypothesis test for normality was 
performed again. The results in Table 16 show that severed additional sites 
are now determined to have normally distributed errors as shown in Figure 
39. Results are given only at the sites which were in the direction of 
motion of the meteorological event propagation since accurate wind speed 
prediction is only accomplished at these sites. Note that the normality 
hypothesis was rejected at 11/14 sites before the storm related variation 
was eliminated from the records and rejected at only 8/14 when the storm 
related variation was eliminated. These results indicate that the variaton 
not directly associated with the storm in a storm front was very often 
normally distributed even if the specific storm induced variation was not 
normally distributed. 
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Table 14 . Table of prediction errors of 2 and 10 minute filtered 
3-10 P.M., DATA SET OF. MAY 2 , 1979 . 


SITE NO. 

ERROR WITH 2 MINUTE 

ERROR WITH 10 MINUTE 


FILTERED DATA SET (MPH) 

FILTERED DATA SET (ITO 

J. 

3.88 

3.24 

2 

3.96 

3.J? 

3 

2.38 

1.67 

4 

5.37 

4.15 

5 

3.42 

2.34 

6 

4.05 

3.01 

8 

2.81 

1.71 

10 

4.63 

3.93 

12 

5.29 

4.32 

13 

4.18 

3.16 

14 

4.22 

2.37 

15 

3.74 

2.52 

16 

4.22 

3.22 

17 

5.70 

4.63 

18 

4.45 

3.72 

19 

3.64 

2.92 

20 

3.27 

2.54 

21 

5.02 

4.43 

22 

5.13 

3.69 

23 

2.79 

2.01 

24 

5-54 

4.63 

25 

5.17 

4.06 

26 

4.94 

3.30 

27 

2.36 

3.43 
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Table 15. Results of hypothesis test for the prediction error at 

2 AND 10 MINUTE FILTERED DATA OF 3-10 F^. OF MAY 4, 1979, 
with Sites 7 1 9, and 11 as references. 


SJIF-NOj _ IQ MINUTE FILTERED 2 HINUT F FTI TFRF[i 


1 

accept 

accept 

2 

accept 

reject 

3 

accept 

reject 

4 

reject 

reject 

5 

accept 

reject 

6 

accept 

accept 

3 

reject 

reject 

10 

reject 

reject 

12 

reject 

reject 

13 

reject 

accept 

14 

accept 

reject 

15 

reject 

reject 

16 

reject 

accept 

17 

reject 

reject 

13 

reject 

reject 

19 

reject 

reject 

20 

accept 

accept 

21 

reject 

reject 

22 

reject 

reject 

23 

reject 

reject 

r\ i , 

^4 

reject 

reject 

25 

reject 

reject 

26 

reject 

accept 

27 

reject 

accept 
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Table 16 . Results of hypothesis for the prediction error at 2 and 
10 MINUTE FILTERED DATA FOR MAY 2 , 1979 , WHEN THE PERIOD 
OF THE STORM ( 7:00 - 8 : 30 ) IS DELETED, 


SITE NO. 

10 MINUTE FILTERED 

2 MINUTE FILTERED 

1 

reject 

accept 

2 

accept 

reject 

3 

accept 

reject 

4 

reject 

accept 

. 5 

accept 

reject 

6 

accept 

reject 

8 

reject 

reject 

10 

reject 

reject 

12 

reject 

reject 

13 

reject 

reject 

14 

accept 

accept 

15 

accept 

accept 

16 

reject 

reject 

17 

reject 

reject 

13 

reject 

reject (accept) 

19 

reject 

reject 

20 

accept 

accept 

21 

reject 

reject 

22 

reject 

accept 

23 

reject 

reject 

24 

reject 

reject 

25 

reject 

accept 

26 

reject 

reject 

27 

reject 

reject 
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Figure 39. Map comparing where the test of hypothesis is accepted 
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The wind speed t° w error was so large when the storm was present that the 
predicted wind power lower limit based on the estimated wind speed lower limit 
would be zero. The wind turbine array would be predicted to be shutdown 
during almost the entire interval that storm front was present using the lower 
limit £(t) - 3 The wind power based on the upper limit of predicted wind 
speed W(t) + 3 c w would indicate that the array power would be at the total 
wind array generation capacity. The actual simulated power for such a record 
on an individual wind turbine would continually cycle between rated capacity 
and shutdown# as shown in Figure 40 for sites 2 and 7. The wind power from an 
array of 81 wind turbines at sites 2 and 7 is shown in Figure 41. The power 
out of the 81 wind turbine array would also cycle from rated array capacity to 
very small levels of power several times as the storm passed through the 
array. Some utilities may attempt to keep the wind turbine array operating 
during such storm induced variations. However# the utility must keep the 
capacity of the array in spinning reserve during the entire period of the 
storm since spinning reserve is based on the lower limit W(t) - 3 a^. Thus# 
the utility would not anticipate power from the array at any time in the 
scheduling of units using the modified unit commitment procedures discussed in 
Section 2. The cyclic power variation out of the array during periods of 
storm induced variation would be compensated by fast responding units under 
feedforward control. The cyclic power variation that could not be handled by 
the response capability of these fast responding units would be clipped by the 
closed loop coordinated blade pitch control or the wind turbines in the array. 
The coordinated blade pitch control would not allow the wind generation level 
out of the array to significantly exceed the magnitude of the generation 
change out of the fast responding units under feedforward control so that 
steam generating units and the system AGC would not have to respond to these 
cyclic storm induced wind power variations. The elimination of the need to 
compensate for wind power fluctuation using regulating units under control 
would increase fuel economy of these units# reduce maintenance on these 
units# and reduce loss of unit life on these units. 

A second method of estimating wind power prediction error directly rather 
than estimating wind speed prediction error and then computing the wind power 
prediction error is proposed. This method does not need to determine the wind 
speed prediction error statistics and does not require determining the 
correlation between wind turbine sites. Moreover, this method does not 
require that the wind power prediction error be normal or time invariant. 

The wind power prediction error at time minutes ahead is predicted using 
an average of the error in predicting wind power over the previous K minutes 

e(t + N) = CLa X Z | P(t - k) - P(t - k) | ? 0.1 P(t) } (7.1) 

K k=1 


if the error is greater than ten percent of the power actually produced by the 
array at time t. The predicted error £(t + N) is never allowed to become 
smaller than .1 P(t) since the short term average error 
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Figure 41. Power from an 81 wind turbine array 


(7.2) 


l I iP(t - K) - P(t - k) 

K k-1 

between the actual P(t - k) and predicted P(t - k) wind array power may become 
small over certain periods and will not certainly reflect the expected error 
at t + H The constant C is chosen as 1 but can range from 0-3 depending on 
the forecasted wind conditions that are expected to effect the wind array at 
t + N. In some cases, the wind power prediction error is expected to increase 
and C is selected to be greater than one if a storm is anticipated to pass 
through the array. In other cases, the error (7.2) may at present be large and 
is expected to decrease as the front or storm passes and C is selected to be 
less than one. 

This method of predicting wind power prediction error does not assume the 
error is time invariant but is slowly time varying. The error band based on 
this error 

P(t + N) - e(t + N) < P(t + N) £ P(t + N) + e(t + N) (7.3) 

does not require that the error to be normal and can be adjusted using C so 
that the future error need not be maintained at C = 3 as in the first method 
but can be set at values between 0 < C < 3 depending on the expected changes 
in wind power prediction error. This method is an on-line procedure and 
predicts the error N minutes ahead rather than using error computed off line 
using the entire actual and predicted. This second method is thus far 
superior to the first method discussed earlier. 

This second wind power prediction error prediction method was tested on 
the actual and predicted wind power records at sites 1-5 for the period 1-6 
pm. on May 2, 1979. Reference sites 23 and 25 were used to produce the wind 
power predicted at sites 1-5. 

The wind prediction power error predicted 30 minutes in advance based on 
a 15 minute average of the wind power prediction error is shown in Figure 42 
for site 1-5. The band (73) generally contain the actual wind power 
P(t + N). The band for site 1 and site 5 has cyclic variation due to the fact 
that the error P(t) - P(t) has large cyclic variations. Since e(t + N) 
never decreases below 0.1 P(t) the error e (t + N) depend on the average 

K 

| £ |P(t - k) - P(t - k) l> > 0.1 p(t) 

K k-1 

producing large values of e(t + N) that may or may not reflect the actual 
error P(t + N) - P(t + N) that occurs at t + N. 

A K = 30 minute average of the error prior to t rather than a K = 10 
minute average reduces and smoothes the large cyclic variations in e (t + N) as 
shown in Figure 43. The N = 30 minute prediction interval is utilized for the 
results in both Figure 42 and 43. Other prediction intervals were chosen and 
the accuracy of the prediction in terms of period out of the band and the rms 
error during this period out of the band were not affected. Thus, a N = 30 
minute prediction interval based on a K = 30 minute average of previous error 
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TIME IN MINUTES! 1 TO 6 PM) 


appeared to produce the best results in terms or predicting error into the 
future and minimizing the cyclic variation in error observed when the average 
of this error was snail. 
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SECTION 8 


CONCLUSIONS AM) FUTCRE RESEARCH 

The purpose of the research is the development of: 

(1) a modified unit commitment; 

(2) a modified generation control; 

(3) a trend wind power predictor required by both the modified unit 
commitment and generation control procedures developed; 

(4) a wind power error predictor required by both the modified unit 
commitment and generation control procedures 

These four developments permit one to answer the following two questions which 
are to be addressed by the research to be conducted within the federal tfiods 

Goecgy Beseaccb floor 12&fclS&Q. 

(1) what is the magnitude of the capacity credit that can be assigned to 
wind energy produced by large arrays based on methods for setting 
and meeting the load following and operating reserve requirements 
within a utility's unit commitment. This magnitude of the 
capacity credit assigned to wind will determine the breakeven 
point in terms of 30 year levelized cost in $/KWH that wind energy 
technology must achieve to warrant large scale implementation by 
utilities. 

(2) Develop a generation control strategy that minimizes the impact of 
large rapid changes in wind array generation that is caused by 
"rotor synchronization" [1, A- 3. HI-201 of all wind turbines in the 
array for large meteorological event wind speed changes. The 
utility's steam turbine generation is slow responding and cannot 
compensate for these large rapid wind generation changes. Attempt- 
ing to force these units to compensate for these large wind generat- 
ion changes would cause cycling in these units that would expend 
significant fuel, increase maintenance, and possibly reduce unit 
reliability and lifetime. A modified generation control is proposed 
in this research based on an hour ahead prediction of wind power 
change. This modified generation control would utilize these slow 
responding large steam turbine units up to a limit imposed by the 
utility; fast responding diesels, hydros and gas turbines that are 
not presently effectively controlled, and wind turbine array con- 
trol of wind power output as a last resort* 

This modified generation control strategy has been developed to allow 
the utility to determine the level of participation its large steam turbine 
units should have in compensating for large wind generation changes. This 
modified generation control was designed so that quick pickup units provide 
the principal compensation for the large wind change. The modified quarter 
hour updated unit commitment strategy would continually unload these quick 
pickup units and replace them with standby economic, peaking, and regulating 
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units. The purpose of maintaining quick pickup generation in spinning 
reserve is to maintain adequate spinning reserve and lead following margins 
to compensate fee large drops in wind generation. If the wind generation 
increase exceeds the allocated combination of the response capability of 
large steam turbines under automatic generation control and the level of 
quick pickup generation connected and loaded (that could be unloaded and 
disconnected in 15 minutes by the feedforward generation control to compen- 
sate for the wind generation increase) , the closed loop wind turbines array 
control would reduce the wind generation rate of change to the level the PGC 
and feedforward control could handle. Por wind generation decreases, the 
AGC would again compensate for the wind generation change up to the 
capability allocated to tracking wind generation change. The feedforward 
control of quick pickup units would then be capable of connecting and fully 
loading all quick pickup generation within the spinning reserve within 
fifteen minutes to compensate for wind generation decrease. The quarter 
hour unit commitment would schedule their connection and the feedforward 
generation control would set the gain on their governor controls so that 
they would be properly loadecl These quick pickup units, once connected, 
would b» controlled utilizing the area control error signal used for regu- 
lating units under AGC. The participation factor on the quick pickup units 
would be adjusted by the feedforward generation control to obtain the de- 
sired generation change out of these quick pickup units and thereby prevent 
units under AGC from exceeding the allocated response capability assigned to 
compensating for wind generation. If the predicted wind generation decrease 
is greater than the combined response capability of the feedforward genera- 
tion control and the allocated response capability of the units under AGC, 
then the closed loop array control attempts to build up a back off reserve 
on tne wind turbines by clipping wind generation below the level that would 
otherwise be produced given the present wind speed at the particular time. 

This bade off reserve is utilized as a cir- so that when the predicted 
drop in wind generation occurs it is not jtr than could be handled by AGC 
and feedforward control within 15 minutes, ihe development of the back off 
reserve is possible due to the hour ahead prediction interval and the fact 
that the level of wind generation to be clipped is based on the oaxiou® 
predicted wind generation decrease (trend wind power minus the error in the 
trend wind power prediciton) and not just the trend wind power change. 

The following accomplishments of this research project are unique: 

(1) the development of wind Speed pcedictiOO for meteorological events 
and turbulence induced variation. Prediction of wind speed based on 
turbulence alone was performed in [2] , but the magnitude of 
turbulence induced variation is to small compared to meteorological 
event variation that it does not require prediction to assure power 
system reliability and econony; 

(2) the development of a method for predicting the ectQC in the wiod 
POWC predictor; 

(3) the development of wind povec pcedictiOO methods* The assessment of 
different wind power prediction methods, the effect of increasing 
the number of wind power prediction sites in the array, the proper 
siting of these prediction sites, and the proper simulation method 
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for producing array power variation from several wind speed 
prediciton sites were all investigated in the research; 

(4) investigation of a modified unit commitment procedure that would 
greatly increase the capacity credit given to wind generation. 
Without trend wind power prediction and wind power error prediction# 
a utility would not be able to connect or disconnect nonwind 
generation in proportion to predicted wind generation increase or 
decrease respectively. Thus# although one could achieve a capacity 
credit based on LQLP calculations# the operation of the utility 
effectively prevented wind generation from serving any load since no 
nonwind generation capacity is displaced by wind generation. The 
modified unit commitment procedure proposed would increase load 
following and spinning reserve proportional to the magnitude of the 
wind power prediction error. The magnitude of the spinning reserve 
increase at any time# which is proportional to wind power prediction 
error at that time# is the amount of the wind generation that is not 
allowed to be counted as meeting Iced due to the lack of perfect 
prediction of wind power variation. Wind power prediction thus 
permits one to provide capacity credit for wind and improvements in 
wind power prediction accuracy increase the capacity credit given to 
wind array power variation. The research performed in this project 
is the only published research on modified unit commitment methods 
that can utilize wind power prediction to modify the 24 hour unit 
commitment based on predicted wind generation changes. 

(5) the development of a generation control strategy based on the one 
developed in [8] but that utilizes both the trend wind power 
predictor and the wind power prediction error predictor for both 
meteorological event and turbulence induoed variation. The 
generation control utilizes the control philosophy in the priority 
use of automatic generation control# feedforward oontrol# and array 
control but incorporates the effects of predicting meteorological 
events and the effects of wind array power prediction error. The 
generation oontrol strategy proposed would satisfy utility 
reliability requirements while simultaneously assuring economic 
operation. Furthermore# the proposed generation oontrol would limit 
the cycling on large steam units that would increase fuel costs# 
increase forced outages# and possibly reduce unit lifetime. 

The methodology has not been integrated into an individual package. 

Thus# the capabilities and performances of the modified unit commitment and 
modified generation oontrol can not be fully quantified and validated and 
should be evaluated in a future research project. 

These four developments of this research project are now briefly reviewed 
to point out the conclusions reached and contributions obtained. 

A modified unit conmitment procedure is developed that is composed of a: 

(1) 24 hour ahead unit commitment that schedules connection and 

disconnection of non wind and wind generation based on a 24 hour 
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ahead load forecast and a 24 ahead forecast of diurnal wind power 
variation; 

(2) quarter hour updated unit commitment that schedules connection and 
disconnection of quick pickup units and standby economic, peaking, 
and regulating units. This quarter hour updated unit commitment 
depends on a one or more hour ahead prediction of trend wind power 
variaiton as wen as a half hour or more ahead prediction error. 

Methods of setting the operating reserve, spinning reserve, unloadahle 
generation reserve, and load following reserve for the 24 hour ahead unit 
commitment are established The operating reserve is based on a LOLP calcula- 
tion assuming the wind turbine array experiences no meteorological events and 
the forecasted diurnal wind power variation is modeled as a load duration 
curve within the normal procedures used to calculate operating reserve. Spin- 
ning reserve, unloadahle generation reserve, and load following reserve are 
developed based on Equations (1-3) where the wind power prediction is the 
diurnal wind power variation and the wind power prediction error 
Qj(k). QJ(k) are assumed to be small. The 24 hour unit commitment schedule 
is intoned to maximize the economic operation of the utility without consi- 
deration of the impacts of meteorological events that can severly affect 
operating reliability and economy. Since the effects of meteorological events 
and sometimes their occurrence, time of arrival, and time of departure can't 
be accurately predicted 24 hours ahead, the 24 hour ahead schedule neglects 
their effects and schedules the system to maximize economic operation and 
maximize the capacity credit available tc wind generation capacity. 

The quarter hour updated unit commitment would schedule connection of 
quick pickup units and a set of economic, peaking, and regulating units on 
standby that would assure the operating reliability of the utility for the 
large meteorogical event induced wind power variations. This quarter hour 
updated unit commitment would obviously modify the most economical unit 
commitment schedule provided by the 24 hour unit commitment but only enough to 
insure operating reliability. Methods for setting operating reserve for the 
quarter hour unit commitment are discussed but no such method has been 
formally developed at this point. Methods for setting spinning reserve, 
unloadahle generation reserve, and load following reserve for the quarter hour 
updated unit commitment are given that utilize the trend wind power predictor 
and the wind power prediction error predictor to determine Qt,(k+j/4-l) 

and QJ(k+j/4-l). The constraints utilized by the quarter hour unit commitment 
in meeting their spinning reserve, unloadahle generation reserve, and load 
following requirements are given in Equations 4-6. The constraints on minimum 
and maximum generation on a unit, and on minimum startup and shutdown periods 
for units are also discussed 

The quarter hour updated unit commitment attempts to maintain sufficient 
quick pickup generation in spinning reserve and load following reserve to 
respond to sudden wind generation decreases. Quick pickup units are unleaded 
and disconnected by the quarter hour unit commitment and replaced by standby 
economic and regulating units in order to increase spinning and load following 
reserve as well as to decrease fuel costs. Quick pickup generation is con- 
nected and loaded if wind generation decreases occur. 
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A trend wind power predictor was developed in this research project. The 
research in this project showed: 

(1) time filtering wind speeds caused significant distortion of the 
maximum, minimum, and average values in viod speeds PtedicitiQO and 
could introduce significant delays; 

(2) time filtering is not required to determine meteogolcgical event 
propagation direction, the reference groups used to predict wind 
speed at prediction sites in the wind array, or propagation delays 
between referenced and prediction sites. This is a change from the 
wind speed prediction method developed in [20] ; 

(3) the reference measurement sites should encircle the wind turbine 
cluster at a distance of at least 100 miles away from all wind 
turbine clusters. Meteorological events can propagate at speeds 
between 0-100 mph and thus a 100 mile separation allows one or more 
hour ahead trend wind power prediction; 

(4) the reference groups should not contain storm cell induced cyclic 
variation because such variation is site specific and time varying. 
Using reference sites with cyclic storm induced variation prevents 
prediction of the trend changes in wind speed that are associated 
with the storm front and can be predicted; 

(5) the reference groups used for prediction should change when the wind 
shift associated with an incoming front first affects a particular 
cluster of wind turbines. The reference group should change from 
one that is in front of the wind array in the direction of 
propagation of the initial meteorological event trr eference sites 
that are in front of the wind turbine array in the propagation 
direction of the incoming event. 

(6) the use of several wind speed reference sites introduces a spatial 
filtering of wind speed variation associated with a meterological 
event. This spatial filtering associated with the wind speed 
prediction is shown to cause the predicted wind power variation to 
exceed the actual wind power produced by the array by as much as 10- 
20 %; 

(7) several wind prediction sites are required to produce accurate wind 
array power estimates. The error utilizing a single wind speed 
prediction site to simulate a 90 wind turbine array could be as 
large as 100% depending on the prediction site selected within that 
wind turbine array. The error could be reduced to 25% if three 
reference sites are used. The larger the number of prediction sites 
the smaller will be the effect of site specific effects and wind 
speed prediction errors of ary prediction site. If the wind speed 
at each wind turbine is not predicted due to the computational 
burden, then one should select prediction sites so that each 
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prediction site is geographically closest to an equal number of wind 
turbine sites. This method of siting wind prediction sites 
minimizes the site specific effects ard error of any one wind 
prediction site on the total array power prediction; 

(8) the study of five different methods of simulating wind array power 
variations indicates that there can be significant differences 
between the results obtained using different methods. These 
differences are minimized as the number of wind prediction sites 
increases. No one method of simulating wind array power variation 
will be most accurate for all wind conditions since the magnitude of 
the error and site specific variation at a wind prediction site will 
vary with the wind conditions. Since each simulation method mini- 
mizes effects of error at specific sites and accentuate error at 
other sites, no one simulation method can give the most accurate 
estimate of true wind array power variation for all wind conditions; 

(9) the magnitude of the wind array power prediction error depends on 
the magnitude of the storm induced cyclic variation and turbulence 
induced wjnd power variation that can not be predicted using the 
trend wind power predictor. The error in the trend wind power 
predictor due to the spatial filtering in the wind speed predictor 
also contributes to wind array power prediction error. This wind 
power prediction error can be the magnitude of the capacity of the 
wind array during storms since the large cyclic variations can cause 
cycling between zero and rated array capacity. The error is large 
because the cyclic variation can not be predicted using this 
methodology. The wind power prediction error can be kept below 

10% - 25% fOr other wi; d conditions if a sufficient number of wind 
power prediction sites are used to simulate the wind array power 
variation. 

A wind power prediction error predictor was also developed in this re- 
search The wind speed prediction error was shown to be a zero mean and 
normal at sites where wind speed prediction is successfully accomplished. The 
wind power prediction error was shown to be slowly time varying. Thus, a wind 
power error predictor was proposed that averages the absolute error between 
the actual array power output and the predicted array power over a 15 or 30 
minute period and uses this error estimate to predict power 30 minutes ahead. 
This predicted error is not allowed to be less than 10% of the wind array 
power output since even though the error may become very small for a period of 
time it does reflect the error that can be expected to occur at some time in 
the future. This wind power prediction error predictor was thoroughly tested. 
The error band around the predicted array power was shown to effectively band 
the actual wind power variation. 

There "-re three major research tasks that can be clearly identified based 
on the results of this research project: 


(1) investigation of improved wind power prediction method, 

(2) develop simulation and evaluation of the performance of the 
modified unit commitment and generation control, 

(3) a cost/benefit stud/ of the wind prediction/ modified generation 
control-unit ocranitment procedures. 

The present research contract has shown that wind array power prediction 
can be performed with an accuracy fo 10-20% for non storm meteorological event 
conditions but may have errors of 100% of array capacity during storms. The 
wind array power prediction required wind speed prediction at several sites 
within the wind turbine array if the accuracy of the wind array power predic- 
tion error is to be minimized. Utilizing several wind prediction sites would 
increase the computational requirements. The computational requirement would 
increase directly proportional to the number of wind speed prediction sites 
used and inc. eases the number of wind speed measurement sites. An increase in 
the number of reference wind measurement sites that encircle the array and 
increase in the number of prediction sites in the array linearly increases the 
cost of land, cost of towers, cost of sensors, cost of the data acquisition 
system, and the cost of the computer system required to implement the wind 
power prediction. Alternate- methods of acquiring wind speed and direction at 
the reference and prediction sites might be investigated. Use of additional 
meteorological information such as pressure, temperature, radar for storm 
related information may increase the accuracy of the estimates of wind spe*. ^ 
for non-storm meteorological events and permit accurate prediction of storm 
induced wind speed variation This research is needed to further improve 
the accuracy of wind array power prediction for both storm and non-storm 
induced variation and reduce the cost of implementation of the wind power 
predictor. 

The second research task is to develop a simulation program that can 
evaluate and further develop and refine the modified unit commitment and 
generation aontrol strategy. The simulation program could then evaluate the 
performance of the modified unit commitment/ generation control on several 
different types of utilities with different wind conditions as done in the 
General Electric study [8]. Additional utilities would be studied that did 
not carry such large load following and spinning reserve capability. 

A computer program would be developed that would simulate both the modi- 
fied unit commitment and generation control procedures developed in this 
research project. A program for computing the operating reserve for the 
quarter hour up-dated unit commitment operating reserve based on the procedure 
outlined in Section 23 would be developed. The modified unit commitment 
procedure would determine operating reserve, spinning reserve, unloadable 
generation, and load following requirements based on the trend wind power 
prediction and the wind power prediction error predictor. The modified unit 
commitment would include constraints on minimum and maximum generation level; 
minimum start-up and shutdown times; and operating reserve, spinning reserve, 
load following, and unloadable generation constraints. 

The computer program for simulating the generation control would be based 
on the program developed earlier at Michigan State for evaluating the 
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performance of present automatic generation control in handling wind genera- 
tion change. Hie program is similar to the program developed for EPRL, This 
computer program would determine whether flGC should solely handle the wind 
generation changer whether feed forward control of quick pickup generation is 
required/ and finally whether closed loop array control of wind array power is 
required based on predicted trend wind power and the predicted error in trend 
wind power change. The feed forward generation control based on the area 
control error would be developed and tested. The closed loop array control 
for reducing wind power rate of change would be evaluated. 

The objectives of this research task is to determine whether the modified 
unit commitment/generation control can 

(1) maintain operating reliability/ 

(2) minimize the fuel and maintainance costs utilizing unit commitment 
and economic dispatch programs/ 

(3) maintain the impact of wind generation change on system PCC below 
the levels specified as part of the generation control procedure. 

The third task is to perform a cost benefit analysis of the combined wind 
power prediction/modified unit generation control. The cost benefit analysis 
would assess the cost benefit analysis could be performed that would assess 
the cost of the 

(a) land for each met tcwer, 

(b) the tower/ 

(c) the wind speed and direction/ pressure/ temperature/ radar 
measurements/ 

(d) communication link/ 

(e) the computer for calculation of the wind speed prediction for each 
wind prediction site and the simulation of predicted array power 
variation/ 

as a function of the number of met towers measurements or the method used to 
produce an accurate hour ahead predictor. These costs will be compared with 
the performance and the savings produced by 

(a) increased capacity credit/ 

(b) reduced production costs , 

(c) improved operating reliability/ 

(d) reduced regulation costs , 

(e) improved load following control/ 
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(£) reduced maintenance costs on conventional units under automatic 
generation control, 


(g) increased unit life, 

for different wind conditions. This analysis extends the results of Tfesk 1 on 
the error in predicting each wind condition for the trend hour ahead predic- 
tors as a function of the number and location of meteorological towers for 
each by as sess ing the costs for providing these predictors at various accuracy 
levels as a function of numbers of met towers provided for each. Hiis task 
extends the work under task 2 by cumulatively assessing the effectiveness and 
cost of both the unit commitment update and generation control possible with 
the hour ahead predictors for all wind conditions and met tower number and 
location- combinations. A cost versus benefit analysis of each predictor that 
compares its benefits in increased capacity credits and reduced production, 
maintenance and regulation costs against the costs for implementing that 
predictor could be made. This analysis would permit decisions on the impor- 
tance and priority for investment in the hardware for a specific utility. 
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